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I. INTRODUCTION 

The NASA Space Shuttle/Spacelab Electromagnetic Environment 
Experiment (EEE) will measure radiation in the 0.4 to 100 Ghz range for 
the purpose of determining earth-emitted interference levels and RF 
spectrum occupancy. The quantities of data involved will be very large. 
Therefore sophisticated methods of data compression, storage, and 
retrieval must be developed so that the data can be stored efficiently 
and retrieved rapidly in the most useful format for analysis and display. 

Such methods include the use of random access mass storage devices 
coupled with efficient data compression/reconstruction algorithums. 

Under the tasks of. this' contract Techno-Sciences has developed 
and’ implemented a series of computer programs for the analysis, 
compression, storage and retrieval of a data base consisting 
of 28 tape files from the ATS-6 RFIME experiment. There is one mode 2 
file, three mode 3 files and 24 mode 4 files. Two principal, programs 
are now implemented on the NASA Interdata Model 5 computer using a 
Bryant 5 M byte disk as a random access mass storage device. One program 
is designed to input the tape data, output it to disk and plot the data 
points, if desired. The zero order or differential entropy can also be 
calculated. The second program is for data retrieval and anlysis. Arbi- 
trarily chosen frequency averaged segments of any of the files on disk can 
be plotted at arbitrarily chosen frequency increments {in multiples of 
the basic input frequency increment of 10 khz for mode 4 data and 100 -khz 
for mode 2 and 3 data). Data smoothing can also be accomplished if desired to 
remove noise variations by a moving quadratic least squares polynomial fit. 
Finally, RFI emitters can be located and printed out using any of the mode 4 
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files, one of the mode 3 files and either the raw mode 2 file or a 
smoothed mode 2 file generated using the moving average quadratic 
polynomial smoothing program. 

The compression achieved by the present simple fixed format 
implementation is 2.21 over the tape storage format. An overall 
compression of four to one is easily achievable by variable length, 
compact storage methods which can be further extended to 6.4 to 
one using a simple universa.1 coding/decoding algorithm. 

The succeeding sections describe the data, the methods of data 
analysis, reduction, and coding, and the computer programs developed. 
Complete program flow diagrams and source listings appear in the 
appendices. Also appearing in the appendices is an analysis of the 
R-ratio as a signal detector and an analysis of the quantization 
effects on the data. 
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II. THE DATA BASE 

The data used for the .study we re.-provided by 'Westinghouse on 2 digital 
9 track magnetic tapes. A listing of the original tape and file number 
designations appears in Table II. 1. These data are from the ATS-6 
RFIME C band experiment and are formatted' as shown in Table II. 2 and cover 
the 5925 Mhz to 6425 Mhz frequency band. The principal volume of data 
is mode 4 which is recorded in 10 khz frequency increments. Mode 2 
and mode 3 data on the other hand are recorded in 100 khz increments. 

Only modes 2, 3 and 4 were supplied. The (Contract study was* 
concerned with these modes only. 

In addition to header information each file consists of statistics 
measured for each frequency increment based on a linear envelope detector 
output. 39 12-bit quantized samples of the detector output are taken, 
the mean value computed and recorded as a 12-bit Magnitude without sign) ' 
number and stored in a 16 bit (2 byte) word. The sum of the squares 
of the sampled. values is also computed, but not normalized to 39, resulting 
in a 12 + 12 + 6 = 30 bit number which is stored in a 32 bit (4 byte) 

word. For mode 4 data only, the peak value of' the 39 is also recorded as 

/ 

a 12-bit (2 byte) number. Thus 64 bits storage per band are required. 

The data are written on tape in DEC format and need to be converted 
to Interdata (IBM) compatible format. In particular, byte switching must 
be accomplished to get the most signifi cart byte of each word in the right 
place and DEC floating point numbers must be converted to the Interdata (IBM) 
floating point format. This is accomplished by subroutines developed 
especially for this contract. 
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Some difficulty was encountered in the use of the data. The header, 
pulse detection, attitude, and telemetry information do not appear on the 
tapes in some cases. Because of the marginal interest in these quantities 
for this contract effort, these data were only reformatted and stored on 
disk and not used further except for computer printouts. Furthermore, 
numerous "bad points" appear, apparently due to the method of A/D conversion 
used. These are screened out, when required, by "R-ratio" technique used 
by Westinghouse. An analysis of the R-ratio appears in Appendix A. 
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TABLE I 1.1 


ATS -6 RFIME Data Based Used in Study 


Tape # 

File # 

Mode 

68 

201 

2A 

159 

302' 

3A 

159 

403 

4' 

159 

407 

4 

159 

408 

4 

290 

304 

3A 

290 

305 

3A 

290 

401 

4A 

290 

402 

4A 

290 

403 

4A 

290 

404 

4A 

290 

405 . 

4A 

290 

406 

4A 

290 

407 

4A 

290 

408 

4A 

'290 

409 

4A 

290 

410 

4A 

290 

411 

4A 

290 

412 

4A 

290 

413 

4A 

290 

414 . 

4A. 

290 

415 

4A 

290 

416 

4A 

290 

417 

4A 

290 

418 

4A 

290 

419 

4A 

290 

420 

4A 

290 

421 

4A 
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TABLE I I. 2 

Digital Tape Data File Format 
Modes 2 & 3 


Record #1 
Records # 2-41 
Data Record: 

Bytes 0-767 

Bytes 768-1133 


Header 

Data 


128 6-byte Interlaced 
Mean, Sum Squares Values 

Time, Frequency, 
Attitude, Telemetry Data 


Mode 4 


Record #1 
Records #2-392 
Rocord #393 
Data Record: 

Bytes 0-1023 


Header 

Data Record 

Pulse Detection Data 


128 8-byte Interlaced 
Mean {2 byte) Sum Squares 
(4 bytes). Peak (2 byte) 
Values 


Bytes 1024-1133 


Time, Frequency, Attitude, 
Telemetry Data 
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III. DATA FORMATTING FOR STORAGE/ RETRIEVAL 

The data format described- in Section II is not very effective for 
efficient data storage and data utilization. In particular, it seems 
likely that most data useage will be in decibel quantities rather than 
in the statistics that appear orr the digital tape. It is clear that 

3 quantities must be retained to preserve the 3 degrees of freedom in 
the mean, peak, and sum squares statistics. It is also clear that there 
.is a correlation between the three which can be removed by a better 
choice of coordinates. The resulting 3 uncorrelated coordinates can 
then be coded more efficiently. The three quantities chosen for storage 
are the log of the mean, the log of the peak/mean ratio and the log of the 
normalized sum-of-the squares to the mean squared (inverse R-ratio). 

12 

Note that for 12 bit quantization, the original data range is 1 to 2 
(volts) or 0 to 72.2 db for each of the 3 original quantities. On the 
other hand, for the selected statistics, only the log of the mean has 
a 72.2 db range. For the 39 points used at each frequency, the 
logarithms of the R-ratio and the peak/mean ratios have a maximum 
.range of 1 to 39 or 15.9 db for the R-ratio (a power ratio) and 31.8 db 
for the peak/mean ratio (a voltage ratio). Furthermore, the R-ratio is 
primarily useful for data screening in the range of about .5 to 1 - a 3 db 
useable range. Points with an R-ratio less than .5 to .6 or so are. 
rejected as "bad points". 

From the analyses in Appendices A and B, based on noise considerations, 
it is concluded that 6 bits of quantization for the log of the mean and 

4 bits for each of the other 2 quantities is sufficient to achieve a 
level of quantization noise which is negligible compared with the system 
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noise. Thus a reduction in storage requirement from 64 original data 

bits to 14 stored data bits or 4.57 to one is achievable through quantization 

alone. 

The presently implemented scheme uses 8 bits for each of the quantized 
parameters and a range of 0 to 100 db. 8 bits was chosen for ease of data 
manipulation as the Interdata machine is directly addressable in 8 bit 
bytes. Formatting at the 14 bit level can be achieved however. The 100 db 
range was chosen arbitrarily and can be easily changed. At any rate, the 
the present scheme provides a reduction from 64 data bits to 24 or 2.67 
to one. The rms quantization noise is .11 db. 

A Bryant 5 M byte moveable head disk is used for data storage. 

Each disk sector contains 256 bytes. The 128 frequencies in each tape 
record are stored in 3 x 128 = 384 bytes s 1 .5 sectors. 0.5 sector is 
presently used to contain the frequency, time, attitude and telemetry 
values for each record and can, no doubt be compressed significantly. 

Because of the unreliability of the taped data, it is not known what 
reduction can actually be achieved. At any rate, the original 1134 bytes 
is presently stored in 512 bytes, a reduction of 2.21 to one. 

The- above considerations lead to the conclusion that a compression 
in excess of 4 to one .~rcan- - be readily attained-'" through more efficient 
quantization and formatting, alone. 

Entropy and universal coding studies p©rformed"4ound~t'hat a-- 
set of 3 data points can be further reduced to '-0, average of -rro'-bfts 
or 6.4 original bits to one. A universal coding block of 128 points for 
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each coordinate was used for the study. A simple, effective, universal 
coding method which achieves the 10 bit average is implemented' by 

ft 

sending the maximum and minimum value for each of the parameters • 
across each block together with the 128 quantized vectors, where 
the quantization is effected on the reduced range defined by the 
maximum and minimum values. Because of the reduced range, fewer bits 
are needed to describe values within the range, the actual number of 
bits depending on the range within a particular block, some blocks 
being more variable than others. In order to further control the 
range, bad data points are removed by R-ratio screening with replace- 
ment by an interpolated average value. 
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IV. DESCRIPTION OF THE COMPUTER PROGRAM FOR PLOTTING AND STORAGE 

One of the two principal presently implemented programs on the NASA 
GSFC Interdata Model 5 is used for the plotting of data points, 
quantization studies, data formatting, and disk storage. The data is 
input from digital magnetic tape and optionally output to the Bryant disk 
as described in the previous section. A flow diagram of this program, 
labelled RFCAL4 , and source listing appears in Appendix C. Prior to 
program start, the input tape must be positioned to the start of the 
file to be processed. A series of commands are entered from the system 
teletype in response to program prompts. An example of a command sequence 
appears in Figure IV. 1. The underlined quantities are the operator 
responses. The first program query after program start is whether disk 
output is desired. If it is, the program asks whether or not it is 
desired to initialize the disk pack. If it is, the data is output starting 
with sector #100 (sectors 0-99 are reserved for program storage). If 
initialize is not requested, the disk pack is searched for end of file 
which is designated by a zero in the first 2 bytes of the sector immediately 
following a data file. A data file is of fixed length using the storage 
format described in the last section - 80 sectors for modes 2 and 3, 782 
sectors for mode 4, so that it is possible to find an empty location by 
a rapid scan. 

After output or no output is determined, whether or not it is desired 
to print the record headers and/or the data points is determined. The 
first record header is alv/ays printed for identification/verification 
purposes. Figure IV. 2 is an example of a small portion of the beginning 
of the plotter output when the printing of headers is not selected, but 
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COMMAND SEQUENCE FOR RFCAL4 


COMMENTS 


START 

TTor output 

Y 

Y TO INITIALIZE 
N 

P FOR PRINT DATA 
N 

P FOR PLOT 
N 

TAPE?, FILE?-' 

290 410 
MODE? 

4 

BITS? 

8 

Z FOR ZERO ORDER ENTROPY 
N 

STOP 

EOO 


Start of program 

Enter Y for output to disk 

For no output, enter N 

Enter Y if this is the first file on 

this disk pack, otherwise, enter N 

Enter P if record headers are to be printed 

Enter P to plot the data 

Enter the tape and file # 

15 format 
Enter Mode 
Format II 

Enter quantization bits desired 

Enter Z foV zero order entropy, for 
differential entropy, enter anything else 
Program stop 
End of job 


Figure IV. 1 RFCAL4 Program Command Sequence* to Output Data to Bryant 
Disk without Plotting (Program prompts are not underlined, 
operator responses are underlined) 



Figure IV. 2. RFCAL4 Sample Plotter Output (Tape #159, File #403) 



MBI 
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plotting is selected. Note that the record attitude and telemetry information 

are meaningless. The plotting is on a Varian 514 plotter. The dark 

horizontal lines are caused by defective transistors in the plotter. Note 

that the log of the mean and the log of the peak/mean ratio are plotted 

) 

in db between 0 and 100 db and the R-ratio is plotted as a number between 
0 and one. Several bad points are seen where the R-ratio becomes veiy small. 

The final program query prior to execution is whether zero order or 
differential entropy is to be calculated. Upon program termination 
the selected quantity is printed out and can be used in data compression 
studies. 
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V. DESCRIPTION OF THE COMPUTER PROGRAM FOR RETRIEVAL, PLOTTING AND 
RFI DECTECTION 

The second of the two principal presently implemented programs on 
the NASA GSFC • Interdata Model 5 is used for data retrieval and processing from 
the Bryant disk using the format described in earlier sections as output to 
disk by program RFCAL4 described in the previous section. A flow diagram 
and source listing of this program, labelled RFCAL5, appear in Appendix D. 

Upon program start, a message is printed on the system teletype 
requesting an operator command (type of processing to be done). There 
are four commands presently implemented - "STOP", "LIST", "PLOT", and 
.“FIND". Upon completion of the latter three commands, control returns 
to the start .point. Thus the purpose of the first command is to stop 
execution to allow- for exit from the program. If a non-existent command 
is entered, control also returns to the start point. 

The second command, "LIST", provides a listing of all the data files 
presently on disk by designated tape and file number. There are presently 
one mode 2 file , three mode 3 files and 24 mode 4 files on disk. In 
addition, there is one smoothed mode 2 file produced by RFCAL5 for RFI 
emitter detection (see below). 

The third command, "PLOT" allows one to plot out a selected data file ' 
from disk on the system plotter. The frequency increment in multiples of 
the input increment of 100 khz for modes 2 and 3 and 10 khz for mode 4 can 
be selected. Frequency averaged values are plotted for increments 
larger than the. respective input frequency increments. Figure V.l 

is an example of the plotted portion of a mode 4 file, plotted at frequency 



Figure V.l. Portion of a Mode 4 Plot, 100 khz per point 


3304 ’ 


TAPE = 230 FILE: *= 410 STORY 5ECT0I? = 

3 BITS luf ClUfiNllSflTX&N 
YR = 5 DAY *= 109 TIME = 1432 AND 43.450 SECS 

START FPEO- 5S25.0000 HH2 DELTAF = 0. 1000 '11142, 10.0000 MHZ/GRID LIME 

AT fITUDE VALUES 
-2973. 132.32442 
-45.297335 - 
-3. 434335 


04 O'. 405449 439. 02b244904 13. 410332 -4. 72444339334. 057353 

1 10. 43437092434. 14p78l 120S4. 000000- 1953. 3S009S - 1 10. 013302 
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increments of 100 khz. Because 10 input points go into one plotted point, 
the plotted data appear smoother than the original (compare with figure IV.l) 

The plot routine can also select smoothing. In this case, a 32 point 
least squares quadratic polynomial Is used to smooth the data. The 
polynomial moves 16 points at a time to provide overlap and the 16 best 
fit polynomial values from the center of the interval are selected for plot- 
ting (except the first and last 16 frequency points which obviously must 
come from the end of the interval). The 32 points with 16 points of overlap 
and the choice of a quadratic polynomial were selected on the basis of a 
compromise between accuracy of smoothing and speed of operation. Prior to 
smoothing, bad points are screened out by the R-ratio and replaced by an 
interpolated average. If desired, the smoothed values can be output to a 
new disk file to be used as a smoothed reference for subsequent data 
operations. The mode 2 file in the provided data base has been processed 
in this fashion. A plotted segment of this file before and after smoothing 
appears in figures V.2 and V.3. Note that after smoothing, much of the 
noise is removed. The quantizing noise in the smoothed data is at the 
same level as in the original, but is now obvious because of the sharp 
reduction in the system noise. 

The final input' command is to "HIND". This command requires a mode 4 
file to search for RFI emitters and a mode 2 and a mode 3 file for reference. 
The mode 4 file is scanned sequentially. After every 10th frequency, the 
.next mode 2 and mode 3 frequency is selected. The R-ratio for each point is 
checked against an input threshold and skipped if too small. If the R-ratio 
exceeds the threshold, the excess of the mode 4 detector output over the 
corresponding mode 2 reference is compared with an input, threshold. . If 



Figure V.2. Portion of Mode 2 Data, unsmoothed 
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Figure V.3. Portion of Mode 2 Data, Smoothed by a Quadratic Polynomial 
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the threshold is not exceeded, the program goes to the next point. If 
the threshold is exceeded, an RFI emitter has been detected and the RFI 
level is calculated using the selected mode 3 file as an absolute reference. 
The emitter frequency and level are printed on the system plotter, together 
with the R-ratio and mode 2 and 3 levels. An example of such a printout 
appears in figure V.4. 

Program operation is very fast and is limited in speed primarily by 
the plotter. Hence several thresholds can be selected and tried successive- 
ly. 
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TAPES PROCESS 

EDI 10063/201 

230/ 

304 

230/421 


Mode 3 D5W = 

30. - NOiSE T 

HRnSHO 

LD> 3.00 DB R-RATi.8 

THRESHOLD ; 

FREOc MH2 > 

RFKDB'iO 

R-RA 

T 10 

hfoDE 2*. DB ; 

MTiDfc. 3c DB ) 

5374. 7330 

20. 3463 

0,, 

3141 

40. 6250 

61. 3^81 

6000. 0625 

22. 1002 

0. 

3 1 30 

40. 6250 

6D. 1563 

6004. 5234 

20. 4358 

0. 

7773 

40. 6250 

60. 3375 

6033. 4083 

20. 4303 

0. 

7831 

3 9- 8 4 

58. 203 1 

6033. 5033 

13. 3646 

0. 

7533 

33. 84b 8 

58, 3844 

<£034. G430 

13. 7536 

0. 

7734 

33. S4bS 

58. 3844 

<£034. 1133 

13. 3336 

0 . 

&203 

33. 8438 

58, 5844 

<£034. £03 1 

13. 3363 

u. 

367 2 

33. 8438 

b£. 203 i 

6033. 4 727 

16. 3533 

O. 

7852 

40. 2344 

£2. 1034 

6152. 7617 

21. 5057 

0. 

8828 

41. 7363 

63. 2813 

6 152. 7734 

13. 5374 

0. 

Sb3S 

41. 7363 

63. 2813 

6229. 533 1 

18. 3021 

0. 

7539 

33. 4531 

60. 1560 

6225. 5731 

17. 3330 

0. 

30 86 

33. 4531 

60. 1563 

6225. 64, *4 

13. 7032 

O. 

[7422 

33. 4531 

60, 1563 

6226. 1680 

18. 7163 

0 . 

3203 

33. 4531 

60. 5463 

6226. 62-75 

17. 3562 

G. 

3164 

33. 8438 

60. 5463 

6226. 6332 

18. 6480 

0 . 

[7617 

33. 8438 

60. 5463 

6226. 70 70 

17. 3132 

0. 

3438 

35. 8438 

60, 5463 

6226. 746 i 

13. 2677 

0. 

7813 

33. 8438 

60. 1563 

6226. £28 1 

- 51. 

0. 

7363 

33. 8438 

60, 1563 

6226, fij’jS 

13. 3504 

0. 

7 77 3 

39. 8438 

60. 1563 

6226, 8783 

13. 1303 

0. 

3047 

33. 843 £ 

60, 1563 

6226, Sti6 7 

13. 5533 

0 . 

[7773 

33. 8438 

60* 1563 

6226. 3334 

13. 6573 

0. 

7227 

39. 84b8 

60. 1563 

6226, 3063 

18. 3070 

0. 

7852 

33. 8438 

60. 1563 

62*26. 32d / 

1 8 . 3232 

0* 

7148 

33, 8438 

6u. 1563 

6226, 3375 

13. 0332 

0. 

£875 

33. 8438 

60, 1563 

6226. 3432 

18. 0657 

0. 

[7305 

33. 8438 

60, 5375 

6223, 536 1 

22. 04b i 

0. 

5^79 

33. 45b 1 

60. 5463 

6230. 0073 

26. 3256 

O, 

$805 

33. 4531 

£0. 5465 

6236. 1367 

17. 5656 

0 . 

3164 

33. S4b3 

60. 3375 

6236. 2635 

17. 1336 

0 . 

3242 

33. 84 2 '& 

61. 3281 

6296, 4633 

16. 6413 

0. 

8438 

33. 8438 

61. 7188 

6236, 6172 

17. 6364 

0. 

7363 

33. 8438 

60, 3375 

6236. 7266 

17. 4312 

0. 

7635 

33. 8438 

61. 3281 

6236. 7 46 i 

17. 2 16S 

0. 

3086 

35, 8438 

61. 3281 

6236, S036 

17. 6074 

0. £036 

40. 2344 

61, 3281 

6236. 3375 

17. 9562 

0. 3164 

40. 2344 

6G. 3375 

6237, 0463 ■ 

17. 6705 

O. 

7363 

40. 2344 

61. 3281 

6237. 1787 

17. 3102 

0. 

[7656 

4U. 2344 

61. 3281 

6237. 2461 

17. 3941 

0. 

8125 

40. 2344 

61, 3281 

6237. 3477 

17. 6234 

O. 

3047 

40. 2344 

61. 3281 

6237, 5377 

17. 7522 

0. 

7852 

40, 2344 

61. 3281 

6237. 7363 

18. 353£ 

0. 

13306 

40. 2344 

60. 3375 

6309. 0351 

16. 7344 

0 . ! 

3164 

40. 2344 

62. 1034 

6 J03. 037/ 

17. 1525 

0. !7*0O 

40. 2344 

62. 1094 

6303. 2333 

17. 03*5 

0. 

[7533 

40. 2b**4 

62. 1034 

6303. 277b 

17. 57 3 C 

0. 

8036 

40. 2*44 

6*2, 1034 

6303. 233 1 

17. 16*15 

0. 

f7422 

40. 2b44 

62. 1034 

630S. 2363 

17. 573C 

0. 

30 *6 

40. 2j44 

62. .1034 

6303. 3*38 

16. 435: 

0 . 

8086 

40. 2344 

62. bUOC 

6303. 3477 

16. 2311 

u. 

3359 

40. 2344 

62. 9000 

6:503. 4570 

16. 7225 

0, 

8*231 

40. 23*14 

62. 9000 

6303. 5078 

16. 73 3f 

0 . 

7533 

40. 2b44 

62, 5000 

6 juS. 5*3 1 

16. 1715 

0 . 

3534 

4u. 2344 

62. 9000 

6303. 5566 

17. 632^ 

0. 

8008 

40. 2344 

61. 7188 

6303. 5377 

17. 0360 

0 . 

3281 

40. 2b44 

61.71 88 

6303. 7773 

17. 5273 

0. 

£2Ub 

4U. *2b44 

62. 1034 

6303. 763 i 

17. 24 IS 

0. ' 

3008 

40. 2344 

62. 1034 

6*03. 3261 

16. 8724 

0 . | 

[7963 

40. 2b44 

62. 1034 

6303, 3733 

16. 3512 

O. : 

8008 

40. 2344 

62. 1034 

6310. 0664 

16. 7844 

0 . 

8164 

40. 2344 

62. 1034 

6310. 2363 

16. 3755 

0 . i 

[731 3 

40. 2344 . 

62. 1034 


Figure V.4. Sample RFI Emitter Detection Printout 
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VI.. NEW TECHNOLOGY 

There are no reportable new technology items resulting from the 
work under this contract. The following review activities were per- 
formed to determine any reportable items: 

1. The key technological concepts and ideas studied and implemented 
under the contract effort were identified. These consisted of the 
methods of analysis, coding, and .storage of the EEE data. The extent 

to which these ideas represented new techniques as versus an application 
of known techniques was reviewed. 

2. A review of appropriate published literature to determine 

the uniqueness of the ideas developed under the contract was performed. 

3. Ameeting with the technical officer to discuss the results 
of the contact effort and points (1) and (2) in connection with 
efforts performed at GSFC and under contract with other contractors 
was held. 

As a result of the review activities, it was concluded that there 
were no ideas, discoveries, or improvements or reportable items which 
were first conceived or reduced to practice under the contract. 
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VII. CONCLUSIONS 

The coding,, storage, and retrieval of the electromagnetic 
environment data was studied. A simple method of universal coding 
was found which can provide a 6.4 to 1 reduction in storage over 
the presently used, -fixed format tape methods. Display, storage, 
retrieval, and analysis, computer programs were developed and implemented 
on the GSFC Interdata computer. These methods provide a simple, 
practical and effective way of storing and analysing the EEE data. 
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APPENDIX A 
R-RATIO ANALYSIS 

The R-ratio has been proposed as a method of detecting the presence 
of RFI emitters. .It has been determined empirically that its usefullness 
is questionable in this regard. This can also be shown theoretically 
through elementary probabilistic considerations. Let (R^}be a set of 
39 independent envelope (linear detector) observations. The R-ratio 
is 



If R. is the envelope of a pure Gaussian noise process, then R. is 
Rayleigh distributed: 

P(Rp = % exp(-R? / 2c 2 ) 
o 

where o 2 is the power associated with the in-phase and quadrature 
components of the underlying Gaussian noise process. 

The mean value of is readily seen to be 

E[R.l = Svj2 a. 

p 

Similarly, the average value of is found to be 

14 \ - 2 ° 2 - 

Hence the average R-ratio is approximately 
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For a single CW or angle-modulated signal without noise, the 
R-ratio is obviously = 1. However, for amplitude modulation or multiple 
signal situations, an R-ratio nearer to that for Gaussian noise occurs. 
Consider specifically the simple situation of two equal amplitude angle- 
modulated signals. Then, the average value of R^ is 

e£rJ =-§7 J o cose) 2 +(sin8) 2 J^8 

= 4/ir. 

2 

The average value of R^ is 

e[r?]= 2, 

so that average value of the R-ratio is approximately 



8/ir 2 X .81, 


or almost the same as for pure noise. 

For. 39 points, for Gaussian noise, the standard deviation of the 
R-ratio is determined by the approximate expression 


_L 2 9 f 12 

39 ^^R. - aA/2 )+ a A"/ 2 J 

T9 f - 2 * 2 ) + 2 ° 2 
i=l 

i 39 

£3 ir/4 + 2A/4 -^q* t (R./i /2.o - A/4 ) 
39 i=1 1 



(R?/2ct 2 - 1). 
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From this, one can. readily calculate 




= .03834, 


which is greater than the difference between the average R-ratios for 2 


equal amplitude sinusoidal emitters and noise only. 



APPENDIX B 
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DATA QUANTIZATION ANALYSIS 

As discussed in the main text, the data statistics chosen for 

storage/retrieval were the mean, the R-ratio and the peak/mean ratio. 

Because the ultimate use of the mean and peak values is in decibels, 

these quantities are stored in logarithmic form (log companded) on the 
12 

range 1 to 2 . 

The mean is a 12 bit voltage value. Taking unity as the minimum, 
the maximum range is 0 db to 72.2 db. For 39 points the maximum range 
of the peak to mean ratio is 1 = 0 db to 39 = 31.8 db. The R-ratio 
has the maximum range of 1 to 1/39, which taken as a power quantity is 
0 to -15.0 db. The R-ratio, however, is not used except in the approxi- 
mate range of 1 to .5, smaller values than .5 to .6 or so being used to 
screen out "bad points". Hence the maximum useable range of the R-ratio 
is approximately 0 to -3 db. 

The quantization increment must be chosen so that the quantization 
noise is significantly less than the receiver noise in the data. If a 
data range of some maximum value, M db, exists, and q quantization bits 
are used elementary considerations lead to a quantization one sigma error 
of 



assuming- a uniform quantizer is used. 


From Appendix A, for Gaussian noise, the R-ratio has a mean + standard 



deviation-to-mean ratio of approximately 
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1 + .03834 _ 

a/~ tt/4 


.184 db 


For M = 3 db, for negligible quantization noise, 
3/2 q vHT «.184 


or: 


2 q » 4.71. 

Hence 4 bit or so quantization for the R-ratio suffices. 


A similar analysis can be used for the mean. Following Appendix A, 

=v/tt72 _ <j 
i=l 

39 


i = 




Hence the expected value of the mean plus standard deviation-to-the 
expected value of the mean ratio is 


1 + 


1 - ir/4 
39 (rr/4) 


.698 db 


for an M = 72.2 db range. 


72.2/2 q y 12 >> , 698 
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or: 


2 q » 29.9. 


Thus 6 bit quantization suffices for the mean value. 


The expected value of the peak/mean plus standard deviation-to- 
expected peak/mean ratio can be shown to have a value in excess of 
that of a single point which is 


1 + 


1 - tt/ 4 
ir / 4 


= 1.826 db 


Hence the number of bits assigned to the peak/mean ratio for a 
range of 31 .8 db should satisfy 

31.8/2 q rvT« 1.826 


or 


2 q » 5.03 

Hence 4 bit quantization is sufficient for the peak/mean ratio. 

Thus a total of 6 + 4 + 4 = 14 bits is sufficient for the basic 

quantization of each of the 3 data points at each frequency. 

\ 
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APPENDIX C. 

FLOW CHART AND SOURCE LISTING 
OF 

STORAGE/PLOTTING PROGRAM 
(RFCAL4) 






FLOW DIAGRAM OF RFCAL4 























C -LAB=RFCAL4 

C THIS PROGRAM IS WRITTEN TO INPUT RFIME DATA, 

C DECODE THE FORMATTED I* AT A THRU SUBROUTINE 
C RFIME, PLOT \ HE DATA m THE \ API AN 514 PLOTTER 
C AND QUANTISE THE VALUES TO AN INPUT NUMBER 
C OF BITS. MEAN* DB ). MEAlfi/PEAK RAT I®C WHEN AVAILABLE) 

C R-RATIO ARE QUANTISED. THE P-RATI0 IS BETWEEN 
C 0 AMD 1. THE PEAK/ME AK RATI® AND 
C THE MEAN IS BETWEEN Cft-IOO DB 
C QUANTISED VALUES ARE PLOTTED. MAX BITS IS 
C 15. DIFFERENTIAL ENTRSPY OR 0 ORDER 
C ENTROPY OF THE DATA 15 CALCULATED 
C ON BITS LE 10. 

C DATA CAN BE OUTPUTTED TO DISK 
C FOR THE TAPE^FILE NO. THE MODE NO., 

C GMT, FREQ. ATTITUDE k TELEMETRY 
C DAT A, AND THE DATA POINTS. 

C' DISK HAS LOG REX OF 5$2 BYTES. THUS 
C 112 8 ORIGINAL BYTES G£ INTO 512 
C BYTES. 8 BIT QUANT IS USED FUR OUTPUT. 

C DATA STARTS AT SECTOR 100 SO ONE CYL , 

C IS AVAILABLE FOR THE USUAL FORMAT 
C THE FILES ARE PACKED BY 331 RECORDS 
C FOR MODE 4, 40 RECORDS FOR MODES 28*3. 

C THE START OF THE LAST SECTOR IS HERO TO _ 

C DENOTE END. DATA STARTS FROM THERE UNLESS 
C INITIALISE IS SELECTED- 
C 

IMFLICIT INTEGER** CI-N) 

DIMENSION PDRC2048),PDPC2043) 

DIMENSION I NPC 564 >, AT I T UDC 13 ), PDMC 2048 ) 

DIMENSION ITELEMC 3), IFREQC 2 ), I GMT C 5 ) 

DIMENSION IPLOTC 10), IGDC 10) | 

DI MENSI &N I OUT C £28 > ! 

EQUIVALENCE ( I0UTC t ), NT APE ), < I0UTC 2 ), NFILE >*C IOUTC 3 ), MODE ), 

1 C I OUT C £5 ), I NPC 1>> i , 

EQUIVALENCE C I0UK4 ), IGMTC 1 ) ), C I OUT C 9 ), I FREQC 1 ) )> 

1 C I OUT (11), ATITUDC 1 > >, ( I OUT ( 37 ), ITELENC 1 > ) ^ _ 

DATA IGDC 1 ), IGDC 2>,IGDC 3 ), IGDC 4 ), IGD( 5 )> IGDC £ )/l, 20 1, 401, . 

1 £01,801, 1001/ 

C READ. WRITE. RANDOM FUNCTIONS 
DAT R IREAP, IWRITE/76, 44/ 

DATA PEE, SEE, WHY/I HP, 1H2, 1HY/ 

IMX= 10 1 

TW£>16=2. **16 1 

C SET QUANTISATION PARAMS 
XQUANT -( 2. ** 15 )/100. 

RQUANT -2. 4* 15 
C HERO PROBABILITY ARRAY 
D0 24 1=1,2048 
PDRC I )=0. 

pdpcd-o. 

24 PDMC I ) = 0. 

LASTM^O 
' LASTR=0 
LASTP=0 
WRITEC 0, 34 ) 

34 FORMAT C 15H Y FOR OUTPUT ) 

’ READC 0,6) STORE 
IFC STORE. NE. WHY ) GO TO 103 
C FIND IF INITIALISING 
WRITE< 0,104) 

FORMAT ( 16HY TO INITIALISE ) 

READC 0, 6 > CLEAR 
IRNDM= 100 

IFC CLEAR. EQ. WHY ) - GS TO 103 

CALL SVC 1C I READ, 15, 1ST AT, I0UT, 512, IRNDM) 

IFCI0UTC 1 ). EQ. 0) GS TO 103 
INC=80 

IFC MODE. EQ. 4 )INC=2±33 1 
IRNDM=IRNDM+INC % 

GO TO 105 
CONTINUE 
GET HEADER IN 

CALL RFIMECINP) 

C FIND IF PRINT OF HEADER STUFF IS WANTED 

WRIT EC O, 5 ) f) h 

5 FORMATC 17H P FOR PRINT DATA) 

READC 0, 6 ) PRINT Oh 

6 FORMAT C A 1 ) * 

C FIND IF PLOT OF POINTS IS WANTED 0 O^ 

WRITECO, 30) K <c? 

3D FORMAT C 1 1H P FOR, PLOT ) 

READC 0, 6 ) PLOT IT 
C FIND TAPE. FILE 
WRIT EC 0, 41 ) 

41 FORMATC 16H T APE? , FILE?C IS ) ) 

READC 0, 2 1 )NT APE, NFILE 
C FIND MODE 

WRITEC 0, 1CP 

10 FORMAT C l OH MODE^CIi)) 

READC 0, 1 1 ) MODE 

11 FORMAT C 1611) 

INC = 4 

C SPACING OF DATA DEPENDS ON MODE 

IFC MODE. LE. 3. OR. MODE. EO. 13 ) INC=3 
C PEAKS DDMT - APPEAR IN THESE MODES 
C SET UP PLOTTER 

CALL SETS14( INC- 1, IGD, 6, 0, 0, 100 ) 

CALL PLT 5 1 4C IPL0T , 0, 0, 0 ) 


104 


105 


103 
C 



C SET NUMB OF RECS T0 D0 “ 

LINES«40 

IFC INC, EQ. * )LINES=331 
C GET BIT.5 QUANTISATION IF NO OUTPUT 
IBIT5=8 

IFC STORE. EG. WHY) GO TO 37 
WRITEC 0,22; 

22 FORMAT C 10H BIT5?CI5)) * 

READC 0,21 ) I BITS 
FORMAT C 1GIS) 

CONTINUE 

FIND SHIFT RELATIVE TO 15 BITS 
I QUANT = 2 + 1S-I6IT5) 

C FIND IF SERO ORDER OR DIFF ENTROPY IS WANTED 
WRITEC 0, 31 ) 

31 FORMAT (20H S FOR 0 ORDER EHTPY ) 

READC 0 , 6 ) ENTRY 
m 4 J= 1, LINES 

C GET INPUT RECORD IN ARRAY INP 
CALL REIMEC INP ) 

C SWITCH BYTES IN ARRAY 
CALL 5WITCH< INP ) 

UNPACK TIME, FREQUENCY. ATTITUDE VALUES ^ ^ ^ 

I GMT C 1 )= YR, I GMT C 2 )=DAY, I GMT < 3 )=HR-MIN, I GMT C 4 )=SEC 
I GMT C 5 )= THOUSANDS OF SECS 

IFREQC 1 )=MH2/50« IFREQC 2 )= 10, 000 1 S OF MHS/50 
ATITUDC I ). 1= 1, 13 HAS ATTITUDE PARAMETERS AS PER 
DATA ACQUISITION USER 1 S MANUAL. TELEMETRY VALUES 
MUST 8E MOVED FROM INPUT APRAY T0 ITELEM 
PRINT ON FIRST REC ONLY WHEN PLOT IS WANTED 
DO 35 1=1,3 

ITELEM*: I )=INP( 1 + 548) ' 

CALL FARAMSC IGMT, IFREQ, ATITUD, INP ) 

IFC PRINT. NE. FEE. AND. J. NE. 1 ) GO TO 40 
IFC PLOTIT. EQ. PEE. AND. J. NE. 1 ) GO TO 40 * 

WRITEC 3, 19) MODE, LINES 
FORMAT t 6HGN0DE ,I3,I6,8H RECORDS) 

WRITEC 3, 23) IBITS 

FORMATS 1H , 12, 2 1H BITS OF QUANTISATION) 

WRITEC 3, 7) IGMT 

7 ' FORMAT C 5H YR = ,13, 6H DAY =,I4,7H TIME =, 15, 4H AND, 13, 
1 1H. ,13, 5H SECS y 
WRITEC 3,8) IFREQ 

-8 ' FORMAT < 7H FREQ= , 13, 1H. , 14. 7H MHS/50) 

WRITEC 3, 13) 

VALUES) 


21 

37 

C 


C 

C 

C 

C 

C 

C 

c 

c 

35 


19 

23 


13 


14 

15 


16 

17 


18 

C 128 
40 
C 


FORMAT ( 16H0ATTITUDE 
WRITEC 3, 9) ATITUD 
FORMAT C 1H , 6F12. 6) - 

WRITEC 3, 14) 

FORMATC 17H0TELEMETRY VALUES) ' 

WRIT EC 3, 15) ITELEM 
FORMAT C 1H , 1QI6 ) 

IFC PLOTIT. NE. PEE) GO TO 40 
C SET UP PLOT 

WRITEC 3, 16 * 

FORMAT C 4SH0PL0T OF ME ANC DB ), PEAK/MEAN RATIOC DB ), R-RATI0 ) 
WRITEC 3, 17) 

FORMATC 3H 0. , 17X, 3H20. . 22X, 3H40. , 

1 10X, 3H DB, 9X,3H60. , 22X, 3H3D. ,2 IX, 4H100. ) 

WRITEC 3, 18) 

FORMAT C 3H 0. , 16X, 3H . 2, 22X. 3H . 4, 3X, 7HRATI0 , 6X, 

1 3H . 6* 22X« 3H . 8, 22X, 3H I. ) 

DATA BLOCKS PER TAPE RECORD 
CONTINUE 

HPT POINTS TO OUTPUT STORE LOCATION 
NPT = 1 

C 128 DATA VALUES PER TAPE BLOCK 
DO 1 1111= 1, 128 

C DATA BLOCKS SPACED BY 3 OR 4 DEPENDING ON MODE 
III=INC+CI1II-1 )+l 
C NOW GET LOG OF MEAN 

X=AMAXO( INPC III ), 1 ) 

IFC X. LT . 0, )X=X+TW0 16 

AVG-X 

R=XkX 

X=20. .-fcALOGlOC X ) 

C QUANTISE THE MEAN 
KQ=X*XQURNT + . 5 
KG=KQ/I QUANT 
C SET UP OUTPUT 

INPC NPT )=KQ 
MPT =NPT + 1 
IDF=KQ-LASTM+ 1024 
LASTM=KQ 

IFC ENTRY. EQ. SEE )LA5TM=G 
PDMC IDF )=PI l MC IDF )+ 1. 

C NOW GO BACK FOR PLOT 
X=KQ*IQUANT 
X=X/XQUhNT 
K= 10. *X+ 1. 5 
I PLOT C 1 )=K 

X= AMAX0C INPC III+ 1 ), 1 ) 

IFCX. LT.O. )X=X+TW016 
Y=INPCIII+2) 

IFC Y. LT. 0. )Y= Y+TW016 
X=C TWO 16*X+Y )/33. 

C R RATIO= MEAN SQUARED/SUM SQUARES 1£ 1. 

R=R/X 

C QUANTISE R-RATI0 

I RQ= R * RQU ANT + . S 


ORIGINAL Pi 
OF POOR QUj 
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IRQ=IRQ/IGUANT 
C SET IJP 0UTPUT 

INPCNPT >=IRQ 
NPT-NPT + 1 

IDF=IRQ-LASTR+ 1024 
LASTR-IRQ 

IFC ENTPY. E0. SEE )LA3TR=0 T 
PDRC IDF >=PDRCIDF)+1. 

C SHIFT BACK FSR PLOT 
R=IRQ*IQUANT 
R=R/RQUANT 

C SKIP PEAK IF NfiiT IN THIS MODE 
IFC INC- EQ. 3) G0 T0 12 
X= AMAXOC I MFC 1 1 1 + 3 >- 1.) 

IFC X, LT. 0. >X=X+TW016 
C PKRAT-MEAM/PEAK RATIO 

PKRAT = 20. + AL0G1QC X/AVG ) 

C QUANTISE PEAK PATIO 

I PKRAT = PKRAT* XQUANT+. 5 
IPKRAT=IFKRAT/IGUANT 
C SET UP OUTPUT 

INPCNPT )=IPKRAT 
. IDF=IPKRAT~LASTP+ 1024 
LASTP^IPKRAT 

IFC ENTPY. EQ. SEE) LA3TP=0 
PDPCIDF )=Pr-PCIDF)+U 
C DE-QUANTI2E FOR PL0T 

PKRAT =IPKRAT+IGUANT 
PKRAT = PKRAT /XQUANT 
IPL0TC 3 )=PKRAT * 10. + 1.5 
C INCREMENT OUTPUT COUNTER WHETHER 0R NOT 
C PEAK DATA EXISTS 
12 NPT -NFT + 1 

IR= 1000. *R+ 1. 5 ' 

* IPL0T ( 2 )=IR 

IFC PL0TIT. EQ. PEE ) CALL PLT5 14C IPL0T * INC- 1,1,1) 

1 CONTINUE - 

C PUT OUT VALUES IF DESIRED 

IFC STORE. NE. WHY) G0 TO 4 
C PACK BYTES IN 

CALL PACKC INPC 1 ). 384 ) 

C OUTPUT DAT R 

CALL 5VC1CIWRITE, 15. 1ST AT, I0UT, 512, IRNDM) 

IRNDM=IRNDM+2 
IF CIST AT. EQ. 0) G0 T0 4 
WRITEC 0* 10 i ) 1ST AT 

101 F0RMATC 1GHDISK 5TAT = ,I<5) 

PAUSE 1 

4 CONTINUE 
I OUT C 1 )=0 

: : IFC STORE. EQ. WHY )CALL SVC 1C I WRITE? 15, 1ST AT * I0UT, 512* IRNDM) . 

IFCISTAT. EQ. 0) G0 T0 102 
WRITEC 0, 101 ) 1ST AT 
PAUSE 2 

102 CONTINUE ‘ . 

C GET PR0BS, ENTR&PY 

SM=0. 

SR=0. 

SP=0. 

DS 25 1=1,2048 
SM=SM+PDM(I ) 

SR=SR+PDRC I ) 

25 SP=SP+PDPC I ) - 
HM=G. # 

HR=0. 

HP=0. 

DS 26 1=1,2048 
IFCPDMC I ). EQ. 0. ) G0 TS 27 
PDMCI >= PDMCI )/5M 
HM=HM+PDM( I HALSGC PDMC I ) > 

IFCPDRCI >.'EQ. 0. )GS T0 28 

HRsHR^PDR^I^iALDGC PDRC I ) > 0®’ ^OOR QUALETIS 

IFC PDPC I X EQ. 0. ) GS TS 26 
PDPC I >=PDPC I )/SP 
HP=HP+PDPc I >+ ALSGC PDPC I ) ) 

CONTINUE 
HH=“HM/AL0GC 2. ) 

HR=-HR/AL0G( 2. ) 

HP=-HP/ALGG<2. ) 

IPC ENTRY. EG'. SEE >WFITEC 3, 32) 

IFC ENTPY. NE. SEE )WRITEC 3,33) 

FORMAT C 1SH SER0 ORDER ENTROPY ) 

FORMAT C 22H DIFFERENTIAL ENTROPY ) 

WRITEC 3, 29 > HM, HR, HP 

FORMAT ( 15H MEAN ENTRDPY= .F12. 6. 19H R-RATI0 ENTROPY = , 

1 F12. <S, 1SH PEAK ENTROPY =,F12. 

STOP 


27 


26 


32 

33 

23 


ORIGINAL PAGE IS 


END 

, U OOOOR 
PDR 1 13CR 
PDP 313CR 
INP 51BCR 
ATITUD 5150R 
ITELEM 5134R 
I FREQ 514CR 


I GMT 
I OUT 
NT APE 
NFILE 


5142R 
513CR 
5 1 iCR 
513ER 


EXT FUNC 
REAL VAR 
REAL VAR 
INT2 VAR 
REAL VAR 
I NT 2 VAR 
I NT 2 VAR 
I NT 2 VAR 
TNT2 VAR 
1NT2 VAR 
INT2 VAR 
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APPENDIX D 

FLOW CHART AND SOURCE LISTING ' 

OF 

RETR I E VAL/ PLOTTI NG/RFI DETECTION PROGRAM 


(RFCAL5) 
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START 


ENTER COMMANDED I 


^ CMD = 

X ‘ STOP 1 V 


START SECTOR 
ON DISK=KSTSEC 

, . S , 

PRINT HEADER 
' INFORMATION 


■'TMD^v 

LIST‘D 


^CMD^- 

PLOT‘S- 


_ READ MODE 4 
H RECORD 

|KSTSEC=KSTSEC+2 

~ — zx zzz 

DECODE RECORD, 
DATA=INP4(IIII) 


-^CMD 
X' FIND^> 
^j'tES 

f ENTER MODE 2 
[TAPE#, FILE# 


IIII=1 


1 1 153841 , 


III =1 


START SECTOR 
ON DISK=KST2 

t p — t 


ENTER MODE 3 
TAPE #, FILE # 

^^SK2^; 


START SECTOR 
ON DISK=KST3 


ENTER: 

EIRP REF=EIRP 
(S+N)/N THRESH=THRSH2 
R-RATIO THRESH=RTHRSH 


INITIALIZE COUNTERS, 
i PR08. ARRAYS, ETC. 


ENTER MODE 4 
TAPE #, FILE # 


-"ON 

DISK? 


CO 


READ MODE 2, 
MODE 3 RECORD 


DECODE MODE 2, 
MODE 3 RECORD, 
DATA=INP2f III) ,INP3( III 



FLOW DIAGRAM OF RFCAL5 














37 


-'TODE 4^ 
R-RATIO 
>RTHRSHt 


MODE 4 S+N v T| NO 
^MODE 2 N >THR jH2^> 



1111=1111-^3 [ - 

C4in^4?^ 
JL YES 


KSTSEOSTART SECTOR 

i zzr 

J READ DISK SECTOR 
AT KSTSEC 


PRINT TAPE #, FILE #| 


KSTSEC=START SECTOR OF 
NEXT DATA FILE 


FLOW DIAGRAM OF RFCAL5(continued) 
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ENTER 

TAPE #, FILE # | 

1 

<^oTdiskJ>- 

NO 

’^^T^yes 


START SECTOR 
ON DISK=KSTSEC 




1 OUTPUTS' N'‘ 

- I 


FIND EMPTY FILE 
START SECTOR=IOUTSCl 


OUTPUT^ 1 Y'l 


PRINT HEADER U. 
INFORMATION 


± 


ENTER START 
FREQ, FREQ 
INCREMENT 


KSTSEC=START 


1 


SECTOR OF PLOTTED DATA 

3c. 


READ RECORD) 
FROM DISK. 



FLOW DIAi 
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OF RFCAL5 (continued) 
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C LAB=RFCAL5* 

C THIS PROGRAM IS WRITTEN TO INPU’T RFImE DATA. 

C IN COMPRESSED FORM BY RFCAL4 ON DISK 

C PLOT THE DATA -IN THE VARIftN 514 PLOTTER 

C INPUT DATA IS QUANTISED T Q 8 BITS. ON RECONSTRUCTION. 

C THE R- RATIO IS BETWEEN j 

C*'0 AND 4. THE PEAK /ME AN RATIO ANJ> . 

,C THE MEAN ARE BETWEEN GiklOG DB \ 

C DIFFERENTIAL ENTROPY OR 0 ORDER 
C ENTROPY OF THE DATA IS CALCULATED 
C DISK HAS LOG REC OF 512 BYTES. THUS 
C' 1128 ORIGINAL BYTES GO INTO 512 
C BYTES* i 

C DATA STARTS AT SECTOR 100 SO ONE CYL 
C IS AVAILABLE FOR THE -USUAL FORMAT 
C THE FILES ARE PACKED , BY 331 RECORDS 
C FOR MODE 4a 40 RECORDS FOR MODES * 2^3. 

C THE START OF THE LAST SECTOR IS, EERO TO 
C DENOTE END. * 

C ' . ' 1 ‘ ’ 1 

- .IMPLICIT- INTEGER* 2 CI^N) ( 

DIMENSION PDRC IQ24).PDP( i02 l j4) 

INP4C5£4).ATITUDCil3).PDMC 1024) 

ITELEMC 3 >. IF RECK 2 ). IGMTC 5 ). CMDSC 3) 
IPLOTC 1D).IGDC 10). i ; . * 

INP2<5*4),INP3<5*4> 1 ! ^ * 

I RD4C £23 )> IRD2C 628 )'. IRD3C 623 ) 

INPSAVC So ). INP5C 5.54 ) ! , { 

“ C IRD4C 41), INP5C l l ) ) 1 

< INP2C I ), IRD2C £5 ) >. C INP3< l ). IRD3C 65 ) ) 


DIMENSION 
DIMENSION' 
DIMENSION 
■ DIMENSION 
DIMENSION 
DIMENSION 
EQUIVALENCE 
• EQUIVALENCE 


54 

50 

Si' 


EQUIVALENCE ( IRD4< 1 ). NT APE ). C IRD4C 2 >. NFILE C IRD4C 3 ); MODE 
1 < IRD4C6S).INP4C 1)> ! 

EQUIVALENCE C IRD4C.4 ). I GMT ( 1 > )**CIRD4C 3 ). : IFREQC 1) ). . - x " 

1 C IRD4C 11). ATITUDC 1 J >* CIRD4C 37 >+ ITELEMC* i > > j 

DATA I GDC 1 ). I GDC 2 IGDC 3 ). IGDC 4 ), IGDC 5 ). IGDC £ )/l* 201. 40li 
301. 1001/ .f 

NST ART.-NSMQTH/iGG. 32/ * - . 

CMDSC L). CMDSt 2 ). CMDSC 3p. MAXCMD/4HLIST * 4HPLQT; 4HFINB. 4/ 
CMDSC 4 >/4HST OP/ I | 

IREAD.IWRITE/76, 44/ [ 1 i 

PEE. HEE. WHY/1HP. 1H2. 1HY/ - - ‘ * 

TWO Id* XQUANT. RQUANT. IBttTS, IQUANT/65536, . 327. 6S. 32763. . 8, 
WRIT EC 0* 50 > 

FORMAT C 1 2HC0MMAND? C A4 ) ) 


1 SGI, 
; DATA 
•DATA 
DATA 
DATA 
- DATA 
DATA 


CMD 


128/ 


RE ADC- 0.51 X 
FORMAT C A4 ) * 

DO 52 IGMDp 1. MAXCMD 
IFC CMD. EQ. CMDSC I CMD > ) GO TO 
CONTINUE 
‘WRITEC 0. 55 ) 

. FORMAT C3HN0T FOUND) 

. GO TO 54- 

GO TO c-101. 102. 10X 104 ). I CMD 


53 


52 
'55 

53 

C LIST FILES ON DISK 
101 IRNDM=N5TART 
56 CALL SVC1CIREAD. 15, 1ST AT. I RD4. 6. IRNDM) 
IFC NTAPE.. EG. Q ) GO’ TO 54 L 
, WRITEC 3. 66 > NTAPE. NFILE. IRNDM 
i “ I RNDM=I RNDMf 702 * MAXOC MODE. 3 )-2026 

G& TO 56 1 


I- 


103 
• 57 . 


58 


53' 

60' 


6i 


62 

?\ 


PRINT OUT EMMITERS 


l 


WRITEC 0.57) 

FORMAT c 1 1HM0DE 2 REF? ) 

READC0.21) NTREF2. NFREF2 
IRNDM- NSTART 1 

CALL SVC 1C IREfiD. 15, 1ST AT , IRD4. 6 * IRNDM) • 

IFC NTAPE. EO. 0) GO TO 35 

IFC NTREF2, EQ. NTAPE, AND* NFREF2. EQ, NFILE) GO TO 53 
IRNDM- 1 RNDM+ 702i MAXOC MODE'. 37-2026- . 

GO TO 58 - 
K3T2=IRNDM * 

WRITEC 0. £0 ) ' 

FORMAT C 1 1HM2JDE 3‘ REF?) 

READ c 0. 2 1 )NTREF3. NFREF3. , 

I RNDM= NSTART 

CALL SVCK TREAD. IS, J3TAT..IRD4. 6. IRNDM) 

IFC NTAPE. EG, 0 ) GO TO' 35 

IFC NTREF3. EG, NT APE* AND. NFREF3. EQ. NFILE )G0 TO* 62 
IRNDM-IRNDM+ 702*NAXGC3. MODE,)- 2026 
GO TO £ 1 
K$T3=IRNDH‘ 

WRIT ECO. 7 T7 

FORMATC 3ZHEIRP' REFC DB >- NTHRE5HC DB ). RTHRESH ) 

RE ADC 0 . 43 )EI RP . T HRSH2. RTHRSH 
C EIRP= EIRP OF MODE 3 REF. NTHPESH=THRESHOLD ABOVE NOISE, RTHRSH 
C «. R- RATIO THRESHOLD 

NTHRSH=255. fTHRSH2/100. +. 5 
IRTHSH=255. *RTHRSH+.-5 
C GET QUANTISED EQUIVALENTS 
102 CONTINUE 
C 2ER0 PROBABILITY ARRAY 
DO 24 1=1. 1024 
PDRC I )=0. 

PDPC I )=0. 

24 PDMC I )=0. 





LASTM=0 
■ LASTR=0 
LASTP=0 

6 F0RMATCA1) 

C FIND TAPE, FILE 

33 WRITEC Q, 4 1 ) 

41 FORMAT < 1-5H TAPE?, FILE? C 15 ) ) 

READC 0, 2 1 )NTPIN, NFLIN 
21 F0RMATC ISIS) 

C FIND IT ON DISK 
IRNBM=NSTART 
C. START SECTOR. 

34 CALL SV C 1 C I READ, 1 5- 1 ST AT , IRD4,--£, I RNDM ) 

■ ' IF( NT APE. EQ. 0 ) GO TO 35 i . - 

C. HOT FOUND IF 2ER0 ENCOUNTERED s 

IFC NT APE. EQ. NTPIN. AND. NFLIN, EQ. NFILE) GO TO 36 . 

IBINC=8G i 

. IF< MODE. EQ. 4 )IDINC=7S2 ! 

IRNDH=IRNDM+IDINC . | 

: GO TO .34 . ‘ . . ' 

35 . WRITE? 0, 37 ) : 

37 FORMAT < 19HTAPELFILE NOT FOUND) 

GO TO 38 • . 

3<S . KSTSEC=IRNDM 

IF< ICMD. EQ. 2 )G0 TO 70 
0 BRANCH ON PLOT COMMAND 
C ELSE FIND EMMITERS MODE 
11*10 • 

111=384 ! 

CIII POINTS TO M0DE2£i3 DATAS. INCREMENTS AT 1/10 RATE 
IFRQD=- 25001 ! 

C STARTING FREO INCREMENT 

C PUT OUT HEADER INFO i ...... „ 

• • WRITER 3, 76 )NTREF2, NFREF2, NTREF3, NFREF3, NTPIN, NFLIN 

76 FORMAT ( 1£H IT APES PROCESSED : , 3C 15, 1H/, 13 )> ' 

‘ WRITE(3,7?> EIRP, THRSH2, RTHRSH 

77 FORMAT < 13HM0DE 3 BBW = , F5. 0, 18H, NOISE THRESHOLD, FS. 2, 3H D8 

, 1 20H R- RAT 10 THRESHOLD = , F5. 4 ) 

WRI TEC 3,-78) I ‘ • 

78 FORMAT < 13H FREQC MH2 J, 4X. 8HRFIC DBW), 4X, 7HR-RATI0, 5X, 

1- 1OHM0DE 2( DB ) ,2X, 1OHM0DE 3( DB ) ) 

73- CALL SVCK TREAD, 15, ISTAT,IRD4, 512, KSTSEC) 

KSTSEC=KST3EC+2 * , 

IFC NT APE: NE. NTPIN. 0R. NFILE. NE: NFLIN ) G0 T0 54 
C QUIT 0Nr END - 

CALL UNPACKCINP4,384> 

. D0 74 1111=1,334,3 , 

IFRQD=IFRQD+ 1 
11=11+1 . 

’ ‘ IFC II. LE. 10) <30 T0 72 

'11=0 

.C II .COUNTS, F0R M0DE 2&3D AT A 
' III=III+3 C 

C III P0INT5 Tfi NEXT VALUE 

XFC III. LT. 384 )G0 Tib 72 • « 

. “ CALL 5VCKIREAD, 15, ISTAr, IRD2, 5 12, KST2 ) 

CALL UNPACKC INP2, 384 > I 

C NEW MODE 2S.3DATA * ■ 1 

CALL SVCKIREAD, 15, ISTAT, IRD3, 512, KST3) 

CALL UNPACKC INP3, 384 ) \ .... 

KST 2=RST2+2 

- ' K5T3=KST3+2 

111=1 . - 
72 .CONTINUE 
C <30 THRU H0DE 4 DATA 
* IRQ=INP4C IIII+ l ) 

IFCIRQ. LT.IRTHSH) <20 T0 74 
C“ TEST F0R NOISE 

KQ=INP4C IIII ) ! ' - 

IF< KQ-INP2C III >. LT. NTHRSH )G0 10 74 
C STHER TEST FOR N0I3E 

X=C KQ-INP3C III ) HIQUANT . 

C, RECONSTRUCT EIRP 
‘ X=X/XGUANT - 
' X=X+EIRP 

- * R-FL0ATC IRQ+IQUANT >/ROUANT \ 

R3=FL0AT( INP3C III+ i PIQUANT, )/RQUANT 
-C ADJUST FROM MEAN TO P0UER * i 
X=X+ 10..+ AL0G IOC R3/R ) 

. * . X2=FL0AT ( INP2C III HIQUANT )/|XQUANT 

- X3*= FLOAT < INP3C III HI QUANT )/XQUAN7 OF F 

FREQ=£ 175. +. O1+FL0ATC IFRQD )\ 

URITEC3, 75) FREQ, X, R, X2, X3 , 

75 FORMAT t 1H ,5F12. 4) 

74 . CONTINUE . 

<30 TQ 73 * 

70 IBLK=-9f MAXOC MODE, 3 )+37 | 

C INCREMENTS ARE 100 KHH F&R M0DE! 2,3, 10 KH5 FOR MSDE 4 
XBLK= FLOAT ( I BLK )* . 0 1 
C XBLK=FREQ INCR IN MHS 
WRIT ECO, 42) 

42 FOR MAT C le'HFREQ, DELTAFC MHS) : 

READC 0, 43 ) SFREQ,DELTF 

43 FORMAT ( £FS. 2 ) 

C kF=NUMB OF BLOCKS FORWARD 
C FOR MODE 4 DATA 



KF=*C3FREG-5325. VL ZB 

IFCKF. GE. 0. AND. KF. LT. 391 ) G0 T0 44 

WRITEC 0. 45) 

45 • F0RMATC 12H0UT 0F RANGE ) 

G0 T& 3£ 

44 KF=KF/X6LK 
C-BLfiCK I NCR ADJUSTED F0R M0DE 
IAVG=DELTF^XBLK+. 5 

* IAVG-MAXGCIAVG, 1 ) i 

0 NUMB 0F PM NTS T0 AVG MUST BE A ( T LEAST i 
XAVG=I AVG 

' DELTFR-XAVG*XBLK 
11 F0RMAT C ifel 1 ) . 

C SET UP PLOTTER 

CALL SET514C3,IGD,<£,0,0, 1000 
0 0NE GRID LlHE=i00 PTS- 1GG*DELTFR MH2 *' 

GRIDL- 100, *DELTFR 
CALL PLT514CIPL0T,D,O,Q> 

C SET NUMB 0F REC3 T0 D0 
C K PLACE T0 START 

' IRNDM^KSTSEC+KF+KF i 

C REMEMBER, IBLK- 10 F0R MSiDE 2,3 - IBLK=1 F0R M0DE 4*. 

- LINE3=3£9/TBLK+2 

' LINES^LIN ES-KF i ' . 

C FIND IF HER® 0RDER 0R DIFF ENTR0PY IS WANTED , , 

WRITER 0,31) | 1 i 

-31 ‘ FERMAT < 20 H 2 FOR 0 ORDER SENTRY 

READ( G, £ ) ENTPY 
WRITEC 0* SO ) I 

SO FORMAT C 12HY' FOR SM00TH ) 

RE ADC 0, 6; SMTH J 

OUTPUT “G„ 

■IFC SMTH. NE. WHY. OR, MODE, NE. 2.) G0 T0 32 # 

C FIND IF SMOOTHED TO BE OUTPUTTED 
• WRITEC 0; S3) i 

33 FORMAT C 12HY FOR OUTPUT ) 

. - READ (0,^) OUTPUT .. ‘ 

- IFC 0UTPUT. NE. WHY) G0 T0 32 
C FIND NEW TAPE NUMBER 
' WRITEC 0- 41 ) 

READCO, 21) NEWTP, NEWFL 
C FIND START SECTOR 

IOUTSC-IRNDM ■ * ' , ‘ I 

94 CALL SVC K IREAD, 15, 1ST AT, IfL,<-,.£, I0UTSC ) 

, . IFC IRD2C 1 ). EQ. 0 ) GO T0 92 L 

: ' IOUT SC-I0UTSC+7Q2*MAXQ< 3, IRID2C 3 ) )-202£ 

‘ G0 TO 94. ' * 

92 NPLTAV=rO * * . ‘ • 

C COUNTS TO NUMBER 10 AVERAGE T0Ri PLOTTING , 

MNAVG-0 
l • IRAVG=0 
1 ' IPKPK-0 
-C 2ER0 AVG. VALUES T00 

. - INI TAT = 0 ' ' i 

C ADJUSTS SMOOTHING FOR FIRST TIME THR0UGH. 

C RESET START SECT0R - *, 

m 4 J= 1, LINES 


PAGE 0003 

41 


C GET INPUT "RECORD IN ARRAY INP4 I 


512, IRNDM) 


CALL SVC 1< IREAD, 15, ISTAT^IRI 
IRNDM-IRNDM+2’ 

C SET FOR NEXT ACCESS 

C UNPACK BYTES T0 HW - - . 

CALL UNPACKS IMP 4, 334) 

C IGMT< 1 )-'YR, IGMT< 2)^ DAY, IGMTC 3 >=HR-MIN, IGMTC 4 >=S£C 
C I GMT <5 )= THOUSANDS OF SECS - .• 

C IFREQC 1 )=MH5/50, JFREQC 10, 000’ S ‘0F MH2y50 

C ATIIUDC I -).'!= 1, 13 HAS ATTITUDE PARAMETERS AS PER 
C DATA ACQUISITION USER'S MANUAL.' . 


C PRINT 0N FIRST REC ONLY WHEN PLOT IS WANTED ' - 
IFC J. NE. 1 ) GS T® 40 1 .1 

WRITEC 3, £60 NTAPE, NFILE. KSTSEC 

FORMATC^HTAPE = .I^,7H FILE I£, 15H START SECTOR =,I7> 
WRITEC 3, 23 )IBITS 

FERMAT C 1H , 12, 21H BITS 0F QUANTISATION) 

WRITEC 3, 7) I GMT 

FORMAT < 5H YR = - 13, *H DAY =,I4,7H TIME = ,I5-4H AND,I3, 

L 1H. , 13* 5H SEC3 *r 

FREQ-50. 4< FLEATC IFREQC t ) HFL0ATC IFREQC 2 ) ViGOuG. ) - ' 

.WRITERS) FREQ, DELTFR, GRIDL 

FERMAT C 12HST ART FREQ= . F10. 4, 4H MH2* 10H DELTAF - , 
t F 10. 4, 5H MHS,‘*F 10. 4, 14H k MMS/GRID LINE > 

WRITEC 3, 13) 

FORMAT C l^HOATTITUDE VALUES) 

WRITEC 3, S) ATITUD ‘ , 

FORMAT C 1H **F12. 6) . . . ! 

. WRITEC 3, 14^ ' .* ! 

FORMAT C 1 7H0T ELEMET R Y VALUES!) 

WRITEC 3, 15) ITELEM * ' 

F0RMATC 1H , 101^) 

UP PL0T 
WRITEC 3, 16 ) 

FORMAT C 4EHGPL0T 0F MEANC DB )^ PEAK/MEAN RATI0C D8 ), R-RATI0 ) 
WRITEC 3, 17) i 

FORMAT C 3H O.-, 17X, 3H20. ,22X, 3H40. , 

1 10X, 3H DB, 3X, 3H^0, , 22X, 3HS0. , 21X* 4H100. ) 

WRITEC 3, ia> 


23 


S 


13* 


14~ 

15 

C SET 
1 ^ 

17 
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IS ■ ■ FORMAT* 3H 0. , 1£X,3H .2,22X,3H . 4, 9X, 7HRATI0 ,6X, 
‘ 1 3H . £>22X,3H : S, 22X, 3H i. ); 

C 128 DATA BLOCKS PER TAPE RECORD 
40 CONTINUE 

IFtSMTH. HE. WHX) GS TCI SI 
C' GET AN EXTRA BLOCK FOR OVERLAP ■ 

CALL SVCKIREAD, 15. ISTAT, IRD2, 512, IRNDM) 

CALL UNPACK* INP2, 3S4 ) 

C SMOOTH VALUES 

NSM3=3+NSM0TH/2 
NSM2= NSMOT H/2 
NSM4=NSM2/2 
NSM43=3*NSM4 
DO SO 1= 1, NSM43. 

11=1+384 

INP5< I ;=INPSAV< I ) 

INP4CII )=INP2( I ) 

11=384- NSM43+I 
INPSAV( I )=INP4(II ) • 

C SAVE FOR NEXT- TIME 
30 CONTINUE 

DO S2 1111= 1, 384, NSM3 • 

IF* J. EQ. LINES. AND. Till. GE. 3S4-NSM3 )IHITAT=N5M2 
DO S3 11=1, NSMOTH I 

K=IIII+3*( II+NSM4-INITAT-1 y . . 

IRD2CII )=INPS<K) . | 

IRD3* II >=INPS< K+ 1 ) 

IRD4< II )=INP5< K+2 ) 

IF* INP5* K+ 1 ). GT. 12S) GO-T0 ,83 - 
C FILL IN BAD ONES 

IRD2< II )=< IMPS* K-3 )+INPS( K+3 ))/2 
IRD3* II )=*INP5(K-2)+INP5*K+4>)/2 
IRD4* II >=< INP5* K- 1 )+INP5< K+5 ) )/2 
S3 CONTINUE 

CALL SMOOTH* I RD2, NSMOTH) i ' 

-CALL SMOOTH* I RD3, NSMOTH) 

IF* MODE. EQ. 4 )CALL SMOOTH* IRD4, NSMOTH ! 

DO -84' 11=1', NSM2 , 

C INITAT =0 FIRST TIME,' THEN NSM4 i 

III=II+INITAT 
' K=IIII+3*II-3 ’ 

INP4* KO=IRD2( III ) 

INP4* K+ 1 )=IRD3* III ) 

INP4* K+2 >=IRD4< III ) 

84 ' CONTINUE 

- INITAT=NSM4 

C NOW FIRST BLOCK IS TAKEN CARE .OF 
82 CONTINUE ' ‘ 

IF* OUTPUT..NE. WHX ) GO TO Si . 

C SET UP OUTPUT - . 

CALL SVC l* I READ, 15,-IST AT IRD2, 512, IRNDM- 2 ) 

- ‘ DO 95 1= 1, 384 • . • . , 

C PUT -IN SMOOTH VALUES ' 

95 INP2* I )=INP4< I ) 

C PUT IN TAPE, FILE 
IRD2L1 >=NEWTP 
' IRD2C2 )=NEWFL 
C MUST BE MODE 2 
IRD2* 3 )=2 - 

<*■ V TT 

. CALL PACK<INP2',3S4> 

, CALL SVCK I WRITE* 15# 1ST AT* IRD2* 512* I0UTSC) ’ 

■ I£HJTSC=I&UTSC+2 m , . 

81 m l> 1111=1*384*3 * ‘ 

C QUANTISED MEAN 

KQ=INP4< IIII ) 

* MNAVG=MNAVG+K8 
C UPDATE AVERAGE 

IDF-KQ-LASTM+5 12 
LAST Mi KQ 

IFC ENTPY. EOT* SEE )LASTM=0 
PDMC IDF )=PDMC IDF )+ 1. 

C R RAT 10= MEAN SQUARED/SUM SQUARES LE U 
C QUANTISED R-RATIfc 

IR(?=INP4C IIII+ i X 
IRAVS=IRAVG+IRQ 
C UPDATE AVERAGE 

IDF-IRQ-LASTR+5 12 
LAST R= IRQ 

IFC ENTPY. EQ. SEE )LASTR=Q 
PDRC IDF )=PDRC IDF>+ 1. 

C SKIP PEAK IF N&T IN THIS M0DET 
IFC MOTE. NE. 4 > G0 Ti3 12 
C PKRAT “ PEAK/MEAN RAT IS 

IPKRAT=INP4C IIII +2 > 

IFC IFKRAT. GE. 82 )IPKRAT=0 
C CSRRECTI©M F£R PEAK . LT. MEAN 

IPKPK-MAXQC IPKPK* IPKRAT+KQ ) 

C UPDATE PEAK VALUE 

IDF=IPKRAT-LA3TP+512 
LASTP=IPKRAT 

IFC ENTRY. EQ. SEE > LA3TP=Q 
PDPCIDF >=PDPC IDF )+ 1. 

12 " NPLT A V = NPLT A V + l 

* IFC NPLT AV. LT. IAVG >GSJ TW 1 
NPLTAV-0 
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KQ-FLOAK MNAVG )/XAVG+, 5 
X=*KQ*IQUANT , 

X*X/XQUANT 
IPLOT< i )= 10, *X+ 1, 5 
IRO=FLOAT< IRAVG )/XAVG+. 5 
R~ IRQ* I QUANT 
R=:R/RQUANT 

IPLB>T<2)=1000, *R+ 1* 5 
PKRAT-C IFKPK-KQ )*IQUANT 
PKRAT = PKR AT /XGU ANT 
I PLOT < 3 )=PKRAT* 10. + 1, 5 
MNAVG-0 
IRAVG=0 
‘ IPKPk^O 

- CALL PLT514<XPU3T,3, 1, 1 > 

i continue: 

4 CONTINUE 

C PUT OUT 2ER0 IP 0UTPUTTING SMOOTHED 
NT APE=Q 

IF< OUTPUT. EQ. WHY >CALL SVCKIWRITE, 1S> I ST AT, IRD4, 6 , I0UTSC ) 
C GET PR0BS, ENTROPY * ' 

$M=0. . 

SR=0, 

SP*G. L * 

m 25 1*1, 1024 
SH=3M+PDMCI ) 

SR=SR+PDRC I ) - 
25 SP=3P+PDP<I) 

HH= 0. 

HR=0. 

HP=0. 

DO 26 1*1, 1024 
IF< PDMC I ). EG, 0. ) GO TO 27 
. PDM<I.)*PDMCI )/SM 
HM“HH+PDM( I )*AL0G< PDM< X ) ) 

IFCPDRC I ). EQ.O. )G0 TO 2B 
PDR(I)=PDR<I )/<5R 
HR=HR+PDRC X >*AL«GC PDR< I ) ) 

IFCPDP(X). EG. 0. ) G0 T0 2£ 

F'DP< X )*PDP<I )/5P 
HP=HP^PDP< I )*AL0G( PDPC I ) ) 

CONTINUE 
HM=~HM/AL0G<2. ) 

HR=-HR/AL0G< 2. > 

HP=-HP/AL0G< 2# ) 

IF< ENTRY, EG, SEE )WRITE< 3, 32 X 
IF< ENTRY. NE. SEE >WRITE< 3, 33 % 

FORMAT ( 19H 2ER0 ORDER ENTROPY) ‘ 

F0RMATC 22H DIFFERENTIAL ENTROPY ) 

WRITE<3,23> HM, HR, HP 

FORMAT C 15H MEAN ENTR0PY= ,P12. 13H R- RAT 10 ENTROPY * , 

1 F 12. 6* 15H PEAK ENTROPY *;Fii2. 6) 

GO TO 54 


27 


2 & 


26 


32 

33 

29 


104 


STOP 
- END 
. U OOOOR 

PDR - ICIER 
. PDP‘ 2C1ER 

INP4 3C9ER 
ATITUD 3C32R 
ITELEM oC^dSR 
I FREQ 3C2ER 


I GMT 

INP5 

NT APE 

NFILE 

MODE 

PDM 

CMD3" 

IPL0T' 

IGD 

INP2 

IRD2 

INP3 

IRDT 


3C24R 
- 3C1ER 
3C£ER 
3C1ER 
3C20R 
3C22R 
41G4R 
0020R 
5 10£R 
- 0008R 
51SAR 
51 1AR 
56SZR 
S*02R 


INPSAV 5P.EAR 
NSTART 00 ICR 
NSM0TH 001ER 
MAXCMD G040R 
IREAD 0042R 
I WRITE 0044R 
PEE . 004*R 

SEE 004AR 
WHY 004ER 
TW0U 0052R 
XQUANT OOSoR 
RQUANT 005AR 
IBITS 005ER 

IGUANT OOSGR 
54 G0*2R 

50 007AR 
<*H ’ OOOOR 

51 OOBOR 
CMD 5BAAR 

52 ODDER 


EXT FUNC 
REAL VAR 
REAL VAR 
INT2 VAR 
REAL VAR 
I NT 2 VAR 
I NT 2 VAR 
I NT 2 VAR 
I NT 2 VAR 
I NT 2 VAR 
INT2 VAR 
,INT2 VAR 
I NT 2 VAR 
REAL VAR 
REAL VAR 
I NT 2 VAR 
I NT 2 VAR 
I NT 2 VAR 
XNT2 VAR 
I NT 2 VAR 
I NT 2 VAR 
INT2 VAR 
I NT 2 VAR 
I NT 2 VAR 
-INT2 VAR 
INI 2 VAR 
INT2 VAR 
PEAL VAR 
REAL VAR 
REAL VAR 
REAL VAR 
REAL VAR 
REAL VAR 
I NT 2 VAR 
I NT 2 VAR 
LABEL 
LABEL 
EXT FUNC 
LABEL 
REAL VAR 
LABEL 
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ICMD 
53 

55 
101 
102 
103 
104> 
IRNDM 

56 

SVC1 

ISTAT 

66 

MAXO 

57 
21 

NTREF2 

NFREF2 

58 

35 ' 

53 

K3T2 

60 

NTREF3 

NFREF3 

61 

>£2 

KST3' 

71 
43 
EIRP 
THR3H2 
RTHRSH 
NTHRSH 
. ¥ 

IRTH5H 

24 

I 

LASTM 

LA3TR' 

LASTP. 

6 

38 

41 

NTPIN .. 
NFLIN ‘ 
34 
Z6 

IDINC* 
37 * 

KSTSEC 
70 
' IX 
III - 
IFRQD 

76 

77 

78 

73 

UNPACK 

74 

IIII 

72 
IRQ 
KQ 
X 

. U 
R 

FL0AT 

R3‘ 

AL0G1Q 
X2 * 

X3 

FREQ* 

75 

IBLIC 
XBLK * 

42 

SFREQ 
'J>ELTF 
'KF ' ’ 

44 

45 

IAVG 

XAVG- 

DELTFR 

11 

SETS 14 
GRIDL. 
PUTS 14 
LINES 

31 

EHTP* 

SO 

SHTH 

OUTPUT 

32 


5BAER 
011ER 
0108R 
013£R 
03E4R 
01B2R 
1C14R 
5BB0R 
013ER 
OOOUR 
5BB£R 
UDEGR 
OOOOR 
0 ICAR 
04B4R 
5BCCR 
5BCER 
020ER 
0538R 
0272R 
5BD0R 
0232R 
5BD2R 
5BD4R 
02D*R 
033AR 
5BD4R 
035AR 
OASOR 
5BD8R 
5BDCR 
5BE0R 
5BE4R 
OOOOR' 
5BF2R 
0414R 
5BF4R 
5BFAR 
5BFCR 
5BFER 
0452R 
045AR. 
0472R 
5C00R 
5C02R 
Q4C£R ; 
"0570R 
5C04R 
0550R 
5C12R 
09C£R 
5C 18R ► 
5C1ER 
5C24R 
05EER 
0646R 
.06B6R 
070AR 
OOOOR 
03B0R 
5C36R 
Q602R 
5C3CR, 
5C3ER 
5C40R 
OOOOR 
; 5C44R 
OOOOR 
5CS0R 
OOOOR- 
5C58R 
5CSCR 
5C£0R 
09A0R 
SC6CR 
SC76R 
OAQER 
5C7AR 
5C7ER 
5C82R 
0AC4R 
OAAAR 
5C90R 
5C32R 
5C3*R 
0B1SR 
OOOOR 
5C3ER 
OOOOR 
5CA2R 
069£R 
5CA8R 
OBECR 
5CACR 
5CB0R 
0D4£R 


INT2^m 

LABEL 

LABEL 

LABEL 

LABEL 

LABEL 

LABEL 

INTafyA* 

LABEL 

EXTiCUNC 

INT2RVA£ 

LABEL 

EXT* FUKT 

LABEL 

int|Ws* 

IMT2i^A^ 
LABEL. 
LABEL 
LABEL 
I NT 

LABEL ' 

INT2%*& 

INT2^££ 

LABEL . 

LABEL 

INT2WT 

LABEL 

LABEL rt 

REALv?8S 

REAU^gg. 

REALvW: 

I NT 2 ‘Ylfe 
EXT FUgT 
I NT 2 
LABEL _ 

I NT 2 
I NT 2 
I NT 2 Y# 

. i nt 2 m 

LABEL . - 
LABEL. / 
LABEL _ 
INT2 Y®" 
INT2 
. LABEL 
LABEL _ 

I NT 2 Y®' 

* LABEL 
- INT2 . 
LABEL ^ 

I NT 2. Ysg, 
INT2 YSg 
I NT 2 
LABEL 
LABEL - 
LABEL ' 
LABEL ' 
EXT -FLttC, 
LABEL _ 

i nt 2 vse: 

LABEL _ 
INT2 .Y3& 

I NT 2 V3& 
REAL-VSf. 
EXT FWg: 
REAL 
EXT FLJg. 
REAL Yg 
'EXT FL*£~ 
REAL Yg 
REAL 9H6* 
REAL ¥>& 
LABEL _ 
INT2 YSR 
REAL 
LABEL _ 
REAL Y?St 
REAL WR 
INT2 
• LABEL 
LABEL _ 
I NT 2 
REAL Yg 
REAL ¥8R 
LABEL' ^ 
EXT FUHC 
REAL ¥g 
EXT FtfKC 
I NT 2 YSR. 
LABEL 
REAL VAR 
LABEL vrt 
REAL YAg 
. REAL VAR 
LABEL 



PAGE 0007 
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NEWTP 
NEWFL 
IfiUTSC 
94 - 
NPLTAV 
.MNAVG 
IRAVG 
IPKPK 
J MI TAT 

4 
J 

40 

23 

7 

5 

13 
9 

14 

15 ' 
IS ’ 

17 

IS * 

SI 

NS M3 
NS M2 
N3H4 
NSM43 
SO 
82 
S3 
K 

SM&0TH 

84 

95 

PACK 

1 

IDF 
12 - *• ' 

IPKRAT 

PKRAT 

SH 

SR 

SP 

25 

HM 

HR 

HP 

2S 

27 

AL0G 

28 

32 

33 ‘ 
29 

;.S 

.-V' 


0CS2R 
5CB4R 
5CBSR 
5CBSR 
0CF2R 
5C8AR 
5CBCR 
5CBER 
5CC0R 
5CC2R 
1SSSR 
5CC4R 
USSR 
0E34R 
0E7CR 
0F2AR 
0F9SR 
0FD4R 
OFFCR 
103SR 
105ER, 
10ACR 
I10CR 
15SER 
5CCER 
5CD0R 
5CD2R 
5CD4R 
1242R 
14CAR 
13EAR 
5CDSR 
OOOOR 
14B0R 
1514R 
OOOOR 
1S5SR 
5CDCR. 
172AR 
5CE2R' 
5CF0R 
5CF4R 
5CFSR 
5CFCR 
18F2R 
5D00R 
-5D04R 
5D08R 
1A3ER. 
‘ 19B2R 
OOOOR 
1A2SR 
1B50R 
• 1BSCR 
1BBCR 
OOOOR 
OOOOR 


LABEL 
I NT 2 VAR 
I NT 2 VAR 
TNT 2 VAR 
LABEL 
INT2 VAR 
I NT 2 VAR 
I NT 2 VAR 
INT2 VAR 
I NT 2 VAR 
LABEL 
INT2 VAR 
LABEL 
LABEL 
LABEL 
LABEL 
LABEL 
LABEL 
LABEL 
LABEL ‘ 
LABEL 
LABEL 
LABEL 
LABEL 
I NT 2 VAR 
INT2 VAR 
I NT 2 VAR 
I NT 2 VAR 
LABEL, 
LABEL 
LABEL 
I NT 2 VAR 
„ EXT FUNC 
LABEL 
LABEL 
EXT FUNC 
LABEL 
INT2 VAR 
LABEL 
I NT 2 VAR 
REAL- VAR 
REAL VAR 
REAL VAR 
REAL VAR 
LABEL 
REAL- VAR 
REAL VAR 
REAL VAR 
LABEL 
LABEL 
EXT FUNC 
' LABEL 
, LABEL 
LABEL ‘ 
LABEL 
EXT FUNC' 
EXT FUNC 


DOOQ ERRORS 


45 


ORIGINAL PAGE IS 
OE POOR QUALITY 
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PAGE 0001 


$LAB= SMOOTH 

QUADRATIC SMOOTHING PROGRAM 
MAKES BEST MS FIT TO A QUADRATIC 
LDD 1/8/7 6 


SUBROUTINE SMOOT H< IDAT, N) 

C IDAT = INPUT ARRAY, N= NUMBER OF POINTS 
C SMOOTHED VALLES ARE PUT BACK IN IDAT 
IMPLICIT INTEGERS (I-N) 

INTEGER* 2 I DAT , N 
DIMENSION IDAK 1 ) 

XN=N 

C XM= K B AR, XQ= ( K- KBAR )# * 2 BAR 
XM=( XN+ 1. >/2. 

XQ=XM«<XN-i. >/£. 

A=0. 

B=*0. 

C=0. 

DC=0. 

DO 1 K=* 1# N 
X=*IDAI( K ) 

XK=FL0AT < K >- XM 

A=A+X 

B=B+X*XK 

XK=XK*XK-XQ 

DC=DC+XK*XK 

1 C=C+X*XK 
A=A/XN 
B=BAXCHXN) 

C-C/DC 

DO 2 K=1,N 
XK=FL0AT < K >-XM 
- X=A+B*XK+C+< XK*XK-XQ ) 

2 IDAT(K)=X+SIGN(. 5, X) 

RETURN 


END 


SNtyfcTH 

0024R 

FUNC/3UB 

5M0C3TH 

D1D£R 

FUNC VftR 

. Q 

OOOOR 

FXT FUNC 

• P 

OOOOR 

EXT FUNC 

IDAT 

002AR 

FfcRM PAR 

N 

002CR 

F0RM PAR 

*‘XN 

01DER 

REAL* VftR 

. W 

OOOOR 

EXT FUNC 

XM 

01E2R 

REAL VAR 

XQ 

OtEER 

REAL VftR 

ft 

OlFoR 

REAL VftR 

. B ’ 

OlFER 

REAL VAR 

C 

02G2R 

REAL VftR 

DC 

G20£R 

REAL VAR 

1 

0GF4R 

LABEL 

K 

020AR 

I NT Z VftR 


QRIGTNM) page is 
OF POOR QUALITY 



<jW 


PAGE 0005 

X Q20CR REAL VAR ^ 

XK 0210R PEAL VAR 

FL0AT OOOOR EXT FUNC 

Ol&SR LABEL 
IGN OOOOR EXT FUNC 

. ¥ OOOOR EXT FUNC 

0000 ERRORS 
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BSEBSSa PR0GRAM T0 D0 lAb UNDER FORTRAN C0NTR0L LDD 12/75 PAGE 1 

CALL IS * 

* CALL SVC 1C ICMD, LU, 1ST AT * ISTART, IBYTES, IRANDOM) 

% where: 

i ICND= SVCi CMDC BYTE ) 

* LU= L0G UNIT 

* 1ST AT = ST AT US RETURNED 
\ 1ST ART = START ADDR 

4 IBYTES=6YTES T0 XFR 

+ IRAND0M- START SECT0R FDR DIRECT DISK ACCESS 




* LDD 

12/75 





'■ *0PT 

LAB=SVC 1 



□QOOR 


EXTRN 

.0 


OOOQR 


SVCI 

ENTRY 

SVCI 


□QOOR 

DOSO 

0052R 

STM 

3> REG 

SAVE EM 

G004R 

D19F 

0000 


LM 

S,OC IS) 

GET ADDRS 

OOOSR 

279E 


SIS 

9, 14 

CK NUMB ARGS 

OOOAR 

233o 


BSS 

m 


OOOCR 

C8B0 

3333 


LHI 

11,C'33‘ 


001 OR 

4 1F0 
OOOOF 


BAL 

15,. 0 

SEND ERR MESS 

0014R 

E 130 


DC 

X'E130* 

QUIT 

00 1 toR 

48 A A 
0000 

BK 

LH 

10, OC 10) 

GET CMD 

001AR 

D2A0 

00£QR 


STB 

10, SVC 


001ER 

4SBB 

0000 

D2B0 


LH 

11, OC 11 ) 


0022R 


STB 

11,SVC+1 

SET LU 

0Q£1R 

002£R 

40D0 

00S4R 

4ADE 


STH 

13, SVC+4 

SET START ADDR ' 

002AR 


AH 

13, OC 14) , 

* 

0000 

002ER 

27D1 


SIS 

13, i 

END ADDR 

0Q30R 

40D0 

0066R 


STH 

13, SVC+£ 

SET IT 

0Q34R 

48FF 

0000 


LH 

45, OC 15) 

i 

0038R 

40F0 

0QS8R 


ST.H 

15, SVC+S 

SET RAND0M IF ANY 

003CR 

E110 

OO^OR 


SVC 

3, SVC 

DO BPER • 

0040R 

4830 

00£2R 


LH 

3, SVC+2 

GET STAT 

0044R 

409C 
0000 
D 130 


STH 

9, OC 12 ) * 

RETURN STAT 

0048R 


' LM 

’9, REG . 

REST0RE 

0052R 

004CR 

4AFF 

.0000 


AH 

15, OC 15) 


OOSOR 

030F 


BR 

15 

RETURN 

OOS2R 


REG 

’ DS 

14 


00£0R 

0000 

0000 

0000 

0000 

0000 

SVC 

DC * 

0 , 0 , 0 , 0 , 0 


006AR 



END 




t&25!S2S PROGRAM TO D0 IVB UNDER FORTRAN CONTROL LDD 12/75 

m ERRORS 



ORIGINAE PAGE IS 
OB POOR QUALITY 


PAGE 2 



so 


bssesesi 


PROGRAM TO PACK BYTES FROM HV FOR FORTRAN LDD 12T7E 
+ call is: 

* CAi_L PACKCI3TART,IBYTES> 

* where :- 

* 1ST ARTIST ART ADBR 

* IBYTES=NUM6 OF BYTES 


OOQOR 

OOOOR 


*GPT 
4 


LAB^PhCK 

EXTRN . 0 
ENTRY PACK 


OOOOR 

DOCO 

PACK . 

. STM 

12* REG ' ^ 



oorssR ' 





0004R 

D IDF 
0000 


. LM 

13> OC 15) 

. 

OOOSR 

27D6 


- SIS' 

- 13,6 

CK ARGS 

OOOAR 

2336 


B2S 

OK 


OQQCR 

C3B0 

’ * 

LHI 

lUC'ZZ* 

SEND .ERR MESS 


3333 




i 

0010R 

4 IFO 
OOOOF 


. BAL I 

15,. 0 


0014R 

E 130 

• 

DC ! 

- X’ E130 1 

^ QUIT 

0016R 

43FF . 
0000- 

0K 

LH . 

15, OC 15 ) 

BYTES 

QQ1AR 

OAFE 


AHR 

1 15, 14 

END ADDR+1 

00 ICR 

08DE 


LHR : 

! 13, 14 f 

% 

OOIER 

48CD 

0000 

L00P 

LH 

I 12, 0< 13 5 

GET HW ' 

0Q22R 

D2CE 

0000 


STB 

12, OC 14 5 

PACK IT 

G026R 

26E 1 


AIS 

14, 1 

INCR PTRS 

0028R 

26D2 - 


AIS 

* 13,2 

- 

002AR 05EF 


CLHR 

14, 15 

CK D0NE 

002CR 

2087 


BL5 

L00P 


002ER 

D ICO . 
0033R 


LM 

{ 

12, REG 


0032R 

4AFF 

0000 


AH 

* 15, 0< 15 5 

1 


Q036R 

030F 


. 

BR 

- 15 

QUIT. 

0038R 


REG ■ 

DS 

. 8 


0040R 



' END 




NO ERRORS 

♦ .0 

LOOP’ * 
. OK 
* PACK 
REG 


PROGRAM TO PACK BYTES FROM. HW FOR FORTRAN LDD 12775 


00 12R 
OOIER 
0Q16R 
OOOOR 
0038R 


miGw&^¥Em~w 

OF .POOR QUALEESI 


PAGE 


PAGE 2 
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QsaQBsa 


PR0GRAM T0 UNPACK BYTES 10 HU LDD 12/73 

* CALL IS* 

* CALL UNPACK* 1ST ART , -I BYTES > 

* where: 

* 1ST ART = START ADDR 
T IBYTES- NUMB 0r BYTES 

* Cl. E. 1ST ART + 2+ 1 BYTES- 1 * END_ADDR, 


PAGE 


* 

tOPT 

* 

* LDD 

* 


1ST ART* 1=FIRST 
LAB=UNPACK 

12/75 


UNPACKED BYTE L0C 


OOOOR 

OOOOR 



EXTRN 

ENTRY 

• 0 

UNPACK 

SAVE EM 

OOOOR 

DODO 

003AR 

UNPACK 

STM 

13, REG 

0004R 

D IDF 
0000‘ 


; LM 

i 

13, OC 15) 

GET EM 

0008R 

OOOAR 

27D6 

233£ 


SIS 

BSS 

13, £ 
QK 

CK ARCS 

OOOCR 

CSBO 

3333 


LHI 

11,C'33’ 

SEND ERR MESS > 

001 OR 

4 IPO 

oooor 


BAL 

is, .0 ; 

0014R 

E130 


DC 

X* E 130' 

6* QUIT 

00 UR 

4AEF 

0000 

m 

AH 

14, OC 15) 

1 LAST BYTE ADDR+ 1 

. i 

001AR 

Q3DE 


LHR 

13, 14 

1 .. ! 

OOtCR 

4ADF 

0000 


AH 

13, OC 15) 

, LAST HW ADDR+H 

0020R 

27E1 - 

L0E3P 1 

SIS 

14, 1 

DECR BYTE PTR 

0022R 

27D2 


SIS 

13, 2 

DECR HW PTR 

0024R 

D3FE ■ 
0000 


LB 

15, OC 14) 

■ GET A BYTE ‘ 

0028R 

40FD 

0000 


STH 

15* OC 13) 

STSRE'HW 

002CR 

05ED 


CLHR 

14, 13 

CK END ' 

002ER 

QQ30R 

0034R 

2087 

DIDO 

003AR 

4AFF 

0000 


BLS 

LM 

AH 

L00P 
13, REG 

15, OC IS) 

003 SR 

030F 


BR 

15 

RETURN 

003AR 

0040R 


REG 

DS " 
END 

6 t . 

1 


SESEE5E 

m ERRORS 

* . 0 
L00P 
0K 
REG 

* UNPACK 


PREGRAM T0 UNPACK BYTES T0 HW LDD 12/7S 


0012R 
0020R 
00 UR 
Q03AR 
OOOOR 


COOO 


PAGE 2 


!>?' yCX)% °''^ U - 
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' Essuasssaa 

* OOOOR 
-OOGOR 


OOOOR' DG9G 
01F4R 
0004R 43AF' 
OOGO - 
0008R 27A4 
OOGAR. 2337 
OOGCR CSBO \ 
3333 ' 
G010R 4 1FQ 
OOOOF 

0Q14R E 130 
- 0000 
001 j£R 4SAF 
0002 
00 ICR 40A0 
02QSR 
0020R CAAO 
i 04£7 

0Q24R 40AG' 
020AR. 
0023R*E11G 
G2G4R 

QG2CR 4300 
020SR 

0030R 4230 
E flSBMBa 

/.-01P0R 
0034R D 130 
* 0 1F4R 

0G33R 430F . 
0004^ 

003CR DOSG 
► 01F4R. 

0040R 43AF 
; 0000 - 
004 4 R 27A4 
GG44R 2335 - 
Q048R C3hU 
333 -> 
004CR' 44 FO 
0012R 
0050R 48AF 
0002 
0054R 24 B 2 
GQ56R CSC A 
04o* - 
005AR 48DA 
0000^ 
005ER 34 DP 
OOoOR 4GDA 
0000 
0064R Cl AO 
. 005AR 
00£SR D 130 
01F4R 
006CR 430F 
. 0004. 


FACE "1 


PACE 2 
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G070R 

GG74R 

0078R 

G07AR 

007CR 

OQSOR 

0QS4R 

' GOSSR 

OGSCR 

003GR 

OQ94R 

Q03SR 

G03SR 

009CR 


DO 30 

OJF4R 

48AF 

OOOG 

2?AA< 

2337 

CS6G 

3333 

4 1FO 

go4er 
E 130 
GOOD 

0008 

4SFF 

0002'. 

D3AE 

0404 

083ft 

S0A4 

40AF 

oooo 

C48G , 


« 




l 

ARAMS 

STM 

| 3. REG ' 



LH 

10. QMS) 



SIS 

1 10. 10 



BS3 

0K2 



LHI 

' 1 1.0*33* 



BAL 

O 

in 

H 



SVC 

! 3.0 

5 

K2 

LH 

1 14. is5 
\ 


, 

LH 

‘15. 2( 15) 



LB 

!*10. 1023< 


‘ 

LHR * 

; ii. 10 



SRLS 

! 10.4 ■ • 



STH 

: 10.0C 15) 


* 

NHT 

1 h F 1 


• CK ARGS 

ERR MESS 
e»QUIT , 

A< IMP V 
ft< I GMT > 

‘ YR. 100* S DAYS 

t i 

YR'S 

YR STORE 
inn 'A DAYS . 




EA 


D RFIME TAPE, INTO FOR 


RAM S/7 4 


PAGE ■ * 3 


00A0R 

OOA4R 

00A8R 

OGAAR 

00ACR 

OOAER 

00S2R 

OOBoR 

0Q83R 

! 

OOBCR 

GOCOR 
0QC2R 
• 0GC4R 

oocaR 

OOCAR 

GOCCR 

00D0R 

00D2R 

OOD^R 

1J0D8R 

OODAR 

OODER 

GOEOR 

GCE4R* 

OGE£R 

OOESR 

OOECR 

OOEER 

00F2R 

GOF£R 
GOr 3R 

OOF CR 

01 OOP 
G102R 
0104R 
GlOoR 

G1GAR 

G1GCR 


GOOFi 
4CA0 
OiCER 
D3DE 
0405 
03CD * 
3GD4 
OASD 
4CA0 
01CER 
C4C0 
00 GF 
OABC 
40BF * 
0002 
4SDE 
04G2 
08BD 
30BE - 
4CA0 
.OICER 
G3CD 
30CA 
C4C0. ' 
OOOF 
□ABC 
4CA0 
OiCER 
G8CD 
30C7 
C4C0 
0007 
GABO 
4CAG 
OICER 
G8CD 
3GC3 
C4C0 
OuGF 
OABC 
4GBF 
0004 
C-IDO 
0007 
0S8D 
4CA0 
OICER 
4SDE 
0400 

oscd 

90CC 

OABC 

40SF 

0006 

08BD 

SGS3 


MH ; 

i 10. DIO * 

t 

LB 

13. 1023C 

LHR * 
SRLS 
AHR 
MH 

; 1*2. 13’ * 

: 13.4 
- 11. 13 
! 10. DIO 

NHI 

MZ.X'F* 

AHR 

STH 

' 11. 12 
j 11. 2< 15) 

LH 

.! 13. t!026C 

LHR 

SRLS 

MH 

* 11.13 
111. 14 
J -10. D 10 

LHR - 
SRLS 
NHI 

*! 12. 13* 

! 12. 10 
l 12. X 1 F* 

AHR 
MH * 

! U. 12 : 

, 10. DIO . 

LHR ‘ 
SRLS 
NHI 

r i2. 13 : 

1 12. 7 
; 1 2. X 1 7 1 

AHR 
MH - 

. 11. 12 
i 10. DIO 

i * 

LHR ' 
SRLS 
NHI 

*1 12. 13 
I 12.3 

j 

AMR 

STH 

* 11. 12 ■ 
j 1 1.4M5) 

NHI 

I 13. X 1 7 l 

LHR 

MH 

* ! 11. 13 
, 10.D10 • 

LH 

; 13. 1Q24C 

LHR 

SRLS 

AHR 

STH 

4 12. 13 
12. 12 
'11.12 
* 11. 15) 

LHR 

SRLS 

Ml. 13 , 

11. a 


*10 DAY5.1 DAYS 


ACCOM % 

* 10 . 

GET DAY 

i 

ACCOM & . 

ST ©RE. DAY. 

10 HR. 1 HR. 10 MIN. 1 MEN* 10 SEC 

i 

SAVE 
10 HRS • 

* '10 


1 HRS 


ACCOM » 

* 10 

10 MINS 


ACCUM 

8, MAKE R30M FGR NEXT ‘ 
1 MINS 


ACCUM 

STORE- HR- MIN 


10 SECS 


* 10 

SECS, . 1 SECS,. 01 SECS, .001 SECS 


ACCUM SECS*.' 
STORE 

. 1 SECS . 


ORIGINAL PAGE IS- • 
OR POOR QUALM 



* 01G ER 

C112R 

01 MR 
0118R 
GUAR- 
OilCR 

0120R 

oiz^k 

012 6k 


C480 
GOOF 
4CAG * 
01CER 
3301' 
8QC4 
GABO 
4CA0 
G ICIER 
C4P0 
00 Of . 
0 h6D 
*OgF 
0008 


012AR 

012ER 

0132R 

-013£R> 

013AR 

G13CR 

G13ER 

0142R 

0 1 46R 

OMAR 

014CR 

0150R 

0152R 

0 15£R 
0158R 

01 SCR 
015ER 
0 16GR 

0M4R 

OMSR 

GMSR 

016AR 

OMCR 


4SFG 
0200ft 
48FF 
0004 
48DE ' 
040.5 . 
C4Du 
OOFF 
08BD 
S0&4 * 
4CA0 - 
01CER 
CABO *• 
0054 . 
‘C4D0 
OOQF ‘ 
OABD ' 
40SF 
GOQO , 
•2434 
48DE 
0408 
G76B ' 
4CA0 
01CER 
08CD 
91D4 - 
C4CG - 
FOOD- 
* 90CC ; 
OABC 
2731 
2033' * 
40BF 
0002 


* 0170R , C39E"'. 
- “ ‘0414 * 

C174R.4SF0 
020QR 

*0178R 4 8FF 

• ' 0005* 
017CR D3AS 

* * 0001 
0180R C5Au 

0080 ; 

, 0 184R 4863 

. . 0002 , 
Q183R 43C3 
0000 
0 18CR\ C4C0 
' ’ 3F8G. 

01S0R OoDC • 
0132R ,C4D0 * 
0 1 SO - 
0195R 2337 
G19SR CACO 
0080 
01 SCR EAAG 
000 1 

Q1A0R 4300- 
019GR 
G1A4R 90 C 1 
01A5R 05 AC 
G1A8R 4803 
0000 
G1ACR C4C0 
COOG* 
GiBOR 05AC 
G1B2R 40 A F 
0000 
01B5R 40SF 
00G2 

’ 0 IBAR 25F4 
01 BCR 2584 : 
. GIBER C59E 

•ORIGIBMi PAGE IS g - 4S 
OS POOR QUAHm 


FjEf 

D RFIME-TflPE 

* INTO FOR ‘ 


*/ - 

NHI , 

i 1 1. X* F 1 



MH 

; 10, DIO 



LBR ' 

* 12, 13 



, SRLS 

. 12^4 



AHR 

' 1 1. 12* 


[ 

MH 

J 10, *D 10 



, NHI 

1 13> X* F 1 



AHR 

| 11, 13 



STH - 

-I il> 8C 15 > ; 

* 

NOW 

GET FREQ ' - f 



LH . 

j IS* RTS 1 . 


i *, ,, 

LH 

; LS^4( h S') : 



' LH 

.:13,-1030( l 



NHI 

! 13,X-Fp' 

- 



LHR ‘ 

[> i 1- 13 



' SRLS 

VIU 4 

, 


MH 

: lo, d io 



' * * AHI 

1 ii,:ioo 

- 

• 

‘ NHI' 1 

j 13/3Tf 1 


. 

AHR 

|l*lTl3 


• 

STH * 

; 1U0C15) 

■ 


\ LIS 

9“. *4* 



LH ! * 

, 13,.1032C t 



XHR * 

it 

F 

4' 

MH 

jMG, DIO 


' \ '' 

' LHR 

i 12, 13 



SLLS 

rl3,4 



NHI 

^ 12, X'* F00 3 



SRLS 

4 12, 12 


, 

: - AHR ' 

I 11, 12 


- 

- - SIS : 

19, 1 



BN23. 

\ P4* " *. 



* STH 

\ 1 1, 2C 15) 

' -( 

NSW 

GET ATTITUDE 



LHI 

* 3, 1044C 14 


' -* 

LH 

15, R 15 



. LH. 

JE.5C J5.y 

REf 

D RFIME TAPE 

INTO FOR ’ 
i , fc * s , _ 

f 

7 

j 

LB * 

* io, 1C9) 

l 



OHI 

10» a' SO* 



LH 

! fi;ao) 



;lh 

;‘12j CK 9 ) 



NHI 

: 12--X' 3Fi ) 

F 

5 , 

‘ LHR' ^ 

i 13/ 12 



NHI, 

' 13,X‘Q1S) 



B2S 




AHI 

12>X‘3Q‘ 



RRL* 

10> 1 



' B 

1 ps : 1 

F 

5 

SRLS 

! 12^ i 



^HR 

: 10, ia 



LH 

i 12, OC 9 ) 
1 



‘ NHI 

; 12.x* coo-) 



0HR - 

10- 12 


* 

STH 

. 10, 0< 15 ) 


- 

• STH 

i ll.,2<: 15 > 
1 ' - 



ATS 

•15,4 - 



AIS 

(3,4 



CLHI 

} 3, 1G36C 1 ! 


4) 
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PAGE 


rsrA h 


*10 

. .01 SECS 

• 001 SECS 

570RE tOGO 1 S 3 ECS 

REST S3 RE ARG PTR 

ACIFREQ) ' 

FIRST* PART- * . f 

MASK G£<?TD PART - * 

START' ACCOM - 
IQ’S S F MH275Q 
‘ MOTE 0VER 

ALWAYS 100‘S = 1 

j* 1 1 S 0F MH2/5Q 

ACCOM ' ’ .. * 

k STSRE 

SET CNTR 
2D PART ’* 

PEP0 ACCUM 

MOTE 13V ER F0R NEXT r Bci>; CHAR - 
- - » 

* MOTE -0UT 0F LEFT SIDE 


HERE IS BCD CHAR 
ACCUM t 
REDUCE CNT . 

ST0RE, 10,000*5 0F* MHS/50 


GET START ADDR ; 

/ RESTORE ARG-PTR. 

‘ flCAXXJTUPJ <• 

5/74 ' - 

' j » “ 

GET HSP FRAC ' 

SET HIDDEN 1 
' LSP ; . 

1 . *' * 

SET EXP , 

TEST HIDDEN SERIES 

- 2ER0 MEANS DGNE \ 
ADD I TO ZERfi) 

DIVIDE, D&VN FRAC 


PAGE 


EXF/2 FOR M * 3 REP 

OR IN AMS 

NOW GET SIGN ' 


SET AMS 


INCR fiUT PTR 
I NCR IN PTR 
TEST DONE 
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01C2R 4260 
017CR 

01C6R' D 130 
D1F4R 
OiCAR 430F. 
OQCA 

01CER OOOA*' ' 

OIDOR E12C 
01E0R 

01D4R E 120 
01E4R 

01D6R El 20 
* 0202R 
01DCR 4300 
0026R 
OiEGR GOO* 

01 FOR* 
01E4R 0007.. 
immazzm 

oooc 

4S2F 

4F20 

4552 

5220 

0 IFOR 
. 01F4R ‘ 

0200R 
0202R 0001 
02Q4R 4601 
0000 
COOO 
OOOG 

020CR ’ . • 




Iff ERRORS 

. 0 ‘ 

£10 -• 

* G0 * 

INPUT • 
LIST \ 

J L0&P.*-C 
MESS - 
NOG* . 
m 
OKI 
15K2 
P4 
PS 

P 6 ' . 

P7 

PARAMS 
PAUSE * 
R15 
RES 
RFIME 
SWI T CH 
UNPAK 
P4 


0062F: 
01CEF: 
00261’ 
0204F 
01E4F 
- 005AF- 
U1F0F 
r 01D0F: 
r ooiaf 

0Q50F 
0066F: 
,0 15SF 

oisog 

01A4F 
017CF 
0G70F 
G202F 
0200F 
0 1F4F 
OOOOF 
003CF 
OiEOF 
0I56F 


n io 

0G3 


NPGK 
1ST 


ESS 

EG 

15 

tfnUSE 

]|NPUT 


BL 

LM : 
B 

DC - 

SVC 

SVC 

SVC 

B 

DC 
DC 


:P7- 
I 9, REG 

I 

i IOC 15 > 


l 


10 

2* UNPAK . 
2, LIST - 
2, PAUSE 
G0 

6 a MESS ■ 

7 a 12, C 1 1 


RETURN 

UNPACK S* 
PRINT ST AT 
WAIT 

TRT AGAIN. 


RErfD RFIME TAPE! INTO F13R 


140 ERR " 
RAM 5/74. 


DS 
D3 . 
EQU 
DC 
DC 


END 1 
D RFIME TAPE 


■4 > 

\ 14' 

! £-2 * 

1 ' ' ^ 

;x* 480 v*b* 0,0 LU \ RD IN' 


INTO FORTRAN 5/74 




PAGE * 6 

i -* 


PAGE- 1 7 


®RI®1HA£ PVfEJ jg 

ID® jTOOt QUALITY 



BUT PUT RFIME DATA 10/75 

+ OUTPUT 5 >FI HE DAT ft UNDER- FSRTRAN 
I. CALL. CALL IS*, 
i % CALL RFBUT C IGUT ) 

* WHERE : 

¥ IGUT 15 ft A < 43 MW ARRAY* THE 1ST 54 

* HU ARE HEADER IMF®. THE LAST 33* HU ARE I 

i BYTE DATA TO 5E COMPRESSED TO 334 BYTES FOR 
i A TOTAL Or 512 BYTES WHICH IS GUT PUT «5M , 

<r LU 2. THE HEADER IKF13 IS 1 HW FOR TAPE M3.* 
H 1 HU FGR FILE MB- * BM£ HU FOR MODE MB. . I 
v 10 BYTES OF GMT DATA | 

is 4 BYTES OF FREQUENCY DATA ' 

* 52 BYTES G- ATTITUDE DATA 

* SL 13 BYTES OF TELEMETRY DATA 

* 

S LDD 10/7*5 


it LAB-RFGUT 


QGGOR 


' 

ENTRY 

RFGUT 


OOOOR 

. 


EXT RN 

. 0 


OOOOR 

DOSO 

0D7ER 

RF8UT 

STM 

11. REG . 


0004R 

DiBF 


LM, * 

13, OC 15) 

GET ARGS 


0000 




1 

0C03R 

27D4 ‘ 


SIS 

13, 4 

C!< RIGHT NUMBER 

OOOAR 

■233S 


BSS 

m 

GOOCR 

CS30 

3333 


LHI 

1LC33' 

i 

0010R 

41F0 

OOOOr 


BAL 

15* . 0 

SEND ERR MESp 

0014R 

40EO 
00 FAR 

0K 

STH 

14, DATGUT +4' 

SET ADDRS 
FIRST INPUT ku 

001 SR 

CSCE 


LHI 

12, 12SC 14 ) 


0030 




1 

00 ICR 

CAEO 
0 IFF 


A HI 

14,511 

LAST OUTPUT ADDR 
LASl^OUT ADD^ SET 

002GR 

40E0 


STH 

14, DAT0UT+* 


007CR 


0024R 

OSBC 


LHR 

LI, 12 

i 

002£R 

4SDC 

LU 

LH 

13, Of 12 ) 

GET HVf | 

002AR 

0000 

D2DB 


.STB 

13, oc m 

ST0RE fiS BYTE 

OOHER 

0000 
2oB i 


RIS 

n* i 

i 

I NCR BUT ADDR 
IN T0Gi 1 

G030R 

2<?C2 


AIS . 

12,2 

0032R 

05EB 


CLHR 

14,11 

Cl< D0NE 1 

0034R 

GG3£R 

22 S? 
E110 


BNLS. 

SVC 

LUP 

1, DA T SUT . 

SUTPUT 0N LU 2 

0C3AR 

007 ^R 
4SD0 


’ LH 1 

13, DAT OUT +2 

Cl< STAT 

003ER 

0042R 

004SR 

007SR 
4 1F0 
004AR • 
D1B0 
■ 0Q7ER 
430 F 


BAL 
v LM 
B ‘ 

15, CKIT 
U,R£G 
4C 15 > 

RETURN 

G04AR 

0004 

433F 

CKIT 

BE 

OC 15) 

RETURN IF SEE® 

004ER 

0000 

E12C 


SVC* 

2, UNPAK 

ELSE UNPACK STAT 

0052R 

005ER 

E120 


SVC ' 

2, .MESS 

i 

SHEW IT 

OOSoR 

GGo2R' 
El 20 


‘ SVC 

2, PAUSE 

WAIT 

0C5AR 

007 4R 
430F 


B: - 

-12C IS) 

TRY AGAIN' 

0G5ER 

FFF4 
00 Go 

' UNPtfK. 

.DC / 

6* MS 


GGS2R 

007 OR - 
0007 

MESS* 

DC 

7, 14, C 1 1/0, STAT 


0070R 

OOGE 

432F" 

4F20 

5354 

^154 ' 

2020 

- MS *' ' 

DS , 

4 


OOF^r 

0001 ' 

PAUSE 

DC. 

1 * 


C07oR 

3302 

BAT0UT. 

DC 

>:* 3802**0*0* Q 


007ER 

00S8R 

0000 

GOOD 

0000 

REG 

D3 ' 
END 

10^ 
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Ntj EAr.'l'ES 


« V 

00 12R 

C K X 7 

G04mR 

DfsFS UT 

oor^R 

LL'P 

Q026F 

1 1CS3 

00£2R 

MS 

0G70R 

'sK t 

00 i^R 

PAUSE 

00*7 4E 

R£«£ 

0C7ER 

PrStfUT 

00 jOi< 

UNP^K 

OOSER 


GO i 2R 
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